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We review the particle physics aspects of ncutrino-loss double beta decay. This process 
can be mediated by light massive Majorana neutrinos {standard interpretation) or by 
something else {non-standard interpretations). The physics potential of both interpreta- 
tions is summarized and the consequences of future measurements or improved limits on 
the half-life of neutrino-less double beta decay are discussed. We try to cover all proposed 
alternative realizations of the decay, including light sterile neutrinos, supersymmetric or 
left-right symmetric theories, Majorons, and other exotic possibilities. Ways to distin- 
guish the mechanisms from one another are discussed. Experimental and nuclear physics 
aspects are also briefly touched, alternative processes to double beta decay arc discussed, 
and an extensive list of references is provided. 
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1. Introduction: General Aspects of Double Beta Decay and 
Lepton Number Violation 

Neutrino-less double beta decay {01/(3/3) is a process of fundamental importance 
for particle physics. It is defined as the transition of a nucleus into a nucleus with 
proton number larger by two units, and the emission of two electrons^: 

{A,Z) ^ {A,Z + 2) + 2e- (Oz//3/3) . (1) 

There are no leptons in the initial state, but two in the final state. Observation of 
OvPP would therefore show that lepton number, an accidental and classical sym- 
metry of the Standard Model (SM) of particle physics, is violated by Nature. The 
process therefore stands on equal footing with baryon number violation, i.e. proton 
decay. For this reason a huge amount of experimental and theoretical activity is 
pursued in order to detect and predict the process^''^'^'^'^''''^'^'^^'^^'-'^^. 

As well known, the main motivation to search for is the fact that neutrinos 
are, in contrast to the prediction of the SM, massive particles and that basically 
all theories beyond the SM predict them to be Majorana^^ particles. However, as 
we will discuss in this review, there are many other well-motivated particle physics 
scenarios and frameworks that allow for 0i//3(3. Before discussing these aspects, let 
us first give some general comments on lepton number violation. 

Why should we look for Lepton Number Violation (LNV)? The conservation of 
lepton (and baryon) number in the SM is an accidental one at the classical level 
only. In fact, chiral anomalies actually violate this conservation law, and it can 
be shown that the currents associated with baryon and lepton number have non- 
vanishing divergences: 9^ J^'^ = cG^^ G^" ^ 0. Here G^i, is the electroweak gauge 
field strength and = X) 9i 7/^ 9« ; ~ S 7m ^« ■ Though this LNV is not the one 
related to 01/(3(3 or Majorana neutrinos, and the rates of processes associated to it 
are negligible at low temperatures, it shows that lepton number is nothing sacred, 
not even in the Standard Model. Extending the picture from the SM to Grand 
Unified Theories (GUTs), quarks and leptons live together in multiplets, and hence 
both B and L are not expected to be conserved quantities. The combination B~ 
which is conserved in the SM both at the classical and quantum level, often plays an 
important role in GUTs, and is broken at some stage. In the spirit of baryogenesis, 
one needs to require that baryon number is violated, and hence lepton number 
should be violated too. The search for baryon number violation proceeds in proton 
decay, or neutron-anti-neutron oscillation experiments. Lepton number violation is 
investigated in neutrino-less double beta decay experiments, and should be treated 
on the same level as baryon number violation. An observation of LNV would be 
far more fundamental than a "simple" measurement of neutrino properties, which 
is often quoted as the main goal of Oi^/3/3-searchcs. Its implications are far beyond 
that. 

In this review we wish to summarize the particle physics aspects of limits and 
possible measurements of this process. A large number of theories and mechanisms 
to violate lepton number exists, and the most often considered light Majorana 
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Fig. 1. Black-box illustration of neutrino- less double beta decay. 



neutrino exchange (though wen-motivated) is only one possibihty. We should note 
that via the black-box, or Schechter-Vallc, theorem^^ (see also^^), all realizations 
of Eq. (1) are connected to a Majorana neutrino mass. Crossing the process on the 
quark level gives from dd — >■ uue~e~ the relation — )■ udude~e~, and with the 
only input of SU{2)l gauge theory one can couple each ud pair via a. W to the 
electrons, as illustrated in Fig. 1. The result is a Pe^i^e transition, which is nothing 
but a Majorana mass term. Needless to say, this is a tiny mass generated at the 
4-loop level, and too small to explain the neutrino mass scale (or its splitting) as 
observed in oscillation experiments. Naively, one can estimate the contribution to 
neutrino mass as 

(m.)ee ^ ^ 10"'' eV, (2) 



(167r2)4 



where we inserted a factor 1/{16tt'^) for each loop, put for each of the two W 
in the loop, and MeV is the typical mass of the involved electron, up- and down- 
quark. An explicit calculation of the 4-loop diagram with an effective operator as 
the source of yields a very similar number^^. Note that this tiny mass is much 
smaller than the Planck-scale contribution to the Majorana neutrino mass, which 
is u^/Mpi ~ 10~^ eV. There are now two main possibilities: 

(i) the mechanism leading to 0^/3/3 is connected to neutrino oscillation. Here 
there are again two possibilities: 

(ia) there is a direct connection to neutrino oscillation. This is the standard 
mechanism of light neutrino exchange; 

(ib) there is an indirect connection to neutrino oscillation. Examples would 
be heavy neutrino exchange, where the heavy neutrinos arc responsible 
for light neutrino masses via the seesaw mechanism. Another case 
would be -R-parity violating SUSY particles generating Ovfi/S, where 
via loops the same particles generate light neutrino masses; 
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(ii) the mechanism leading to Ovfi^ is not connected to neutrino oscillation. 
The underlying physics of 0v/3f3 could be cither: 

(iia) giving a sub-leading contribution to neutrino mass, maybe i?-parity 
violating SUSY being responsible for only one of the neutrino masses; 

(lib) giving no contribution to neutrino mass, for instance a right-handed 
Higgs triplet in the absence of a Dirac mass matrix for neutrinos. 
Hence only the Schechter-Valle term from Eq. (2) can generate a neu- 
trino mass. 

In both cases (iia) and (lib) we would need another source for neutrino 
mass and oscillation^^. 

As already mentioned, the assumption that massive Majorana neutrinos gener- 
ate 0vf3/3 is presumably the best motivated one, though there are many more. We 
can thus classify the possible interpretations of Oi^/3/3 as follows: 

(1) Standard Interpretation: 

neutrino-less double beta decay is mediated by light and massive Majorana neu- 
trinos (the ones which oscillate) and all other mechanisms potentially leading 
to give negligible or no contribution; 

(2) Non-Standard Interpretations: 

neutrino-less double beta decay is mediated by some other lepton number vio- 
lating physics, and light massive Majorana neutrinos (the ones which oscillate) 
potentially leading to Oiy(3/3 give negligible or no contribution. 

In this review we will consider both cases and aim to discuss all possible real- 
izations of OvPP. In Sections 2 and 3 we will deal with experimental and nuclear 
physics aspects of Oi^/3/3, respectively. The standard interpretation of light neutrino 
exchange is discussed in Section 4, where we summarize in detail our current under- 
standing of neutrino physics and its many aspects which can be tested with Qv^jS. 
Section 5 is devoted to the various non-standard interpretations, such as left-right 
symmetric theories, i?-parity violating supersymmetry, Majorons, and other pro- 
posals. Section 6 deals with possibilities to distinguish the mechanisms from one 
another, and Section 7 is concerned with alternative processes to Oi//3/3. A summary 
is presented in Section 8. For all aspects we provide an extensive list of references. 

2. Experimental aspects 

Neutrino-less double beta decay can only be observed if the usual beta decay is 
energetically forbidden. This is the case for some even-even nuclei (i.e. even proton 
and neutron numbers) , whose ground states are energetically lower than their odd- 
odd neighbors. If the nucleus with atomic number higher by one unit has a smaller 
binding energy (preventing beta decay from occurring), and the nucleus with atomic 
number higher by two units has a larger binding energy, the double beta decay 
process is allowed. In principle 35 nuclei can undergo 0;//3/3, though realistically only 
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nine emerge as interesting candidates and are under investigation in competitive 
experiments, namely ^^Ca, ^'^Ge, ^^Sc, ^^Zr, iooMo, "^cd, "OTc, ^^^Xc, i^^Nd. 
There is no "super isotope", one has to find compromises between the natural 
abundance, reasonably priced enrichment, the association with a well controlled 
experimental technique or the Q-value, because the decay rate for 0^/3/3 goes with 
(except for Majoron emission, see Section 5.5). Table 1 and Fig. 2 give the 
relevant parameters of all 11 isotopes with a Q- value above 2 MeV, including the 
nine most studied isotopes given above. The experimental signal is the sum of 
energy of the two emitted electrons, which should equal the known Q-value. The 
neutrino-less mode has to be distinguished from 2 neutrino double beta decay^^ 

{A,Z) ^ {A,Z + 2) + 2e- + 2D^ (2i^/3^) , (3) 

which experimentally can be an irreducible background for the neutrino-less mode. 
The half-life of 21/(3/3 is typically around 10^^-10^^ yrs (it is important to note 
that the process is allowed in the SM), and has been observed for a number of 
isotopes already, see^^ for a list of results. Obviously, the countless peaks arising 
from natural radioactivity, cosmic ray reactions etc. need to be understood and/or 
the experiments have to be ultrapure and/or heavily shielded. The energy release 
Q should also be large due to the background of natural radioactivity, which drops 
significantly beyond 2.614 MeV, which is the highest significant 7-line in the natural 
decay chains of Uranium and Thorium. In general, the decay rate for 01/(3(3 can be 
factorized as 

r'"' ^GAQ,Z)\MAA,Z)r^^\\ (4) 

where rjx is a function of the particle physics parameters responsible for the decay. 
The nuclear matrix element (NME) Mx{A, Z) depends on the mechanism and the 
nucleus. The term Mx{A,Z)t]x can in fact be a sum of several terms, therefore 

Table 1. Q- value, natural abundance and phase space fac- 
tor G (standard mechanism) for all isotopes with Q > 2 
MeV using ro = 1.2 fm. Values taken from Table 6 of® 
and scaled to = 1.25. Note that there is a misprint in 
Ref.6, which quotes G"" for ^^Omo as 11.3x 10~" yrs'^ 



Isotope 


G [10-1-* yrs-i] 


Q [kcV] 


nat. abund. [%] 


«Ca 


6.35 


4273.7 


0.187 


76 Ge 


0.623 


2039.1 


7.8 


S^Se 


2.70 


2995.5 


9.2 


96 Zr 


5.63 


3347.7 


2.8 


iooMo 


4.36 


3035.0 


9.6 


llOpd 


1.40 


2004.0 


11.8 


116Cd 


4.62 


2809.1 


7.6 


i24Sn 


2.55 


2287.7 


5.6 


i30Te 


4.09 


2530.3 


34.5 


I36xe 


4.31 


2461.9 


8.9 


150Nd 


19.2 


3367.3 


5.6 
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Natural abundance of different O 
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Fig. 2. Natural abundance and phase spaee faetor for all 11 Oi^/3/3-isotopes above Q = 2 MeV. 



including the possibility of destructive or constructive interference, a situation we 
will deal with in Section 6.4. Finally, Gx{Q,Z) is a phase space factor which can 
have dependence on the particle physics. For most of the processes in which only two 
electrons are emitted, the phase space factor can be considered almost independent 
of the mechanism. The biggest effect for G occurs in double beta decay with Majoron 
emission, in which the final state contains one or two additional particles, sec Section 
5.5. Table 2 summarizes the current best limits on the half-life. Ncutrino-lcss double 
beta decay is definitely a rare process. In Table 2 we already quote the limits on 
the effective mass (the particle physics parameter in the standard interpretation) 
from the respective experiments, for which we used a compilation of nuclear matrix 
elements discussed later. 



Table 2. Experimental limits at 90% C.L. on the most interesting isotopes 
for OvPfi. Using the nuclear matrix element ranges from Table 5 we also 
give the maximal and minimal limits on (mee). 



Isotope 


^1/2 


[yrs] 


Experiment 


[eV] 


M'Z. [eV] 


48 Ca 


5.8 X 


10^2 


CANDLES18 


3.55 


9.91 


76 Ge 


1.9 X 


1025 


HDM19 


0.21 


0.53 




1.6 X 


1025 


IGEX20 


0.25 


0.63 


82Se 


3.2 X 


1023 


NEMO-321 


0.85 


2.08 


96 Zr 


9.2 X 


1021 


NEMO-322 


3.97 


14.39 


iooMo 


1.0 X 


1024 


NEMO-321 


0.31 


0.79 


iiscd 


1.7 X 


1023 


SOLOTVIN023 


1.22 


2.30 


i30Te 


2.8 X 


1024 


CUORICIN024 


0.27 


0.57 


I36xe 


5.0 X 


1023 


DAMA25 


0.83 


2.04 


iSONd 


1.8 X 


1022 


NEMO-326 


2.35 


5.08 



Note: The limits on T^J^ from NEMO-3 measurements assume the standard 
light neutrino mechanism. 
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In the past the search rchcd mainly on geo- and radiochemical measurements, 
which are insensitive to the mode of double beta decay, but led to the first obser- 
vation that two neutrino double beta decay occurs in nature^^. Here the approach 
is to identify accumulation of the decay isotope (in particular if it is a noble gas, 
for which mass spectroscopy can be done very precisely) during geological time 
periods in samples which are rich in a double beta decay isotope. In principle, this 
method can be used to test the time-dependence of the parameters associated with 
the mechanism of double beta decay. 

Nowadays only direct methods are applied, based mainly on the observation of 
the two electrons in the form of measuring their total sum energy, which should 
equal the Q-value of the decay. Some experiments have the possibility of tracking 
the individual electrons. There are a number of recent reviews on the experimental 
situation in double beta decay, to which we refer for more details^'^'^^^^^. The 
number of expected events in an experiment can be written as 

N ^\n2aMtNA iT^/2)~'^ , (5) 

where a is the abundance of the isotope, M the used mass, t the time of measure- 
ment and Na is Avogadro's number. The half-life sensitivity depends on whether 
there is background or not^^: 

a M et without background, 

(ri%)-^(x<( rMT (6) 

BAE 



a£ \ ^5— r — with background. 



Here B is the background index with natural units of counts/ (keV kg yr) and AE 
the energy resolution at the peak. In Table 3 we follow the classification proposed by 
A. Guiliani^^ and list some properties of the main up-coming experiments. Table 4 
lists the most developed experiments according to^^. Roughly speaking, at present 



Table 3. Planned experiments categorized according to^^ and the isotope(s) under consideration. 



Name 


Isotope 


high energy 


source = detector 
res. low energy 


calorimetric with 
res. sensit. to event topology 


source ^ detector with 
sensit. to event topology 


CANDLES^s 


«Ca 




/ 






COBRA^o 


"6Cd (and "OTe) 






/ 




CUORE31 




/ 








DCBA^'^ 


isONd 








/ 




I36xe 






/ 




GERDA34 


7606 


/ 








KamLAND-Zen^s 


l36Xe 




/ 






LUCIFER^'S 


**2Sc or ioOMo or "^cd 


/ 








MAJORANA37 


76 Ge 


/ 








MOON^s 


82Se or ioOmo or i^ONd 








/ 


NEXT^s 


l36Xe 






/ 




SNO-l-40 






/ 






SuperNEMO*! 


82 Se or isONd 








/ 


XMASS*^ 


"6xe 




/ 
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Table 4. Sensitivity at 90% C.L. of the seven most developed projects for about three (phase II 
of GERDA and MAJORANA, KamLAND, SNO+) five (EXO, SuperNEMO and CUORE) and 
ten (full-scale GERDA plus MAJORANA) years of measurements. Taken from^^. 



Experiment 


Isotopc 


Mass of 


Sensitivity 


Status 


Start of 






Isotope [kg] 






data-taking 


GERDA 


'^Gc 


18 


3 X 102== 


running 


~ 2011 






40 


2 X 1026 


in progress 


~ 2012 






1000 


6 X 1027 


R&D 


~ 2015 


CUORE 


l30Te 


200 


6.5 X 1026* 


in progress 


~ 2013 








2.1 X 1026" 






MAJORANA 


76 Ge 


30-60 


(1 - 2) X 1026 


in progress 


~ 2013 






1000 


6 X 1027 


R&D 


~ 2015 


EXO 


l36Xe 


200 


6.4 X 1025 


in progress 


~ 2011 






1000 


8 X 1026 


R&D 


~ 2015 


SuperNEMO 


82 Se 


100-200 


(1 - 2) X 1026 


R&D 


~ 2013-2015 


KamLAND- Zen 


"6Xe 


400 


4 X 1026 


in progress 


~ 2011 






1000 


1027 


R&D 


~ 2013-2015 


SNO-f- 


"ONd 


56 


4.5 X 1024 


in progress 


~ 2012 






500 


3 X 1025 


R&D 


~ 2015 



Note: * For a background of 10 'VkcV/kg/yr; ** for a background of 10 2/keV/kg/yr. 



the transition from 10 kg yrs to 100 kg yrs experiments is being made, background 
levels below 10~^ counts/ (keV kg yr) are planned, and half-life sensitivities above 
10^^ yrs arc foreseen. 

The current best values come from the Heidelberg-Moscow^^ experiment, using 
^®Ge enriched Germanium calorimctric detectors. As is well known, part of the col- 
laboration claims observation of the process^'^. at the level of about 2 x 10^^ yrs, 
with a 99.73% C.L. range of (0.7 - 4.2) x 10^^ yrs. This has been criticized by a 
large part of the community'*^, but eventually needs to be tested experimentally. 
In the later part of this review we will discuss limits on lepton number violating 
parameters from Oz^/3/3. As the limit on the half-life of ^^Ge corresponds roughly 
to the claimed signal, one could easily translate the limits of the lepton number 
violating parameters into their values, in case the claim is actually valid. 

A bolomctric experiment, also run in Gran Sasso, was CUORICINO^'*, using 
^•^°Te in the form of Te02 crystals. Similar limits to Heidelberg-Moscow could be 
reached. An experiment with source 7^ detector was NEMO-3'^^, using foils of several 
potential Oj^/3/3-cmitters in a magnetized tracking volume. Here the main point is 
measuring the energy of the individual energies and their angular distribution. This 
approach is of interest in testing different mechanisms for Ov/Sp, as we will discuss 
later. Again, limits of the order of Heidelberg-Moscow were obtained. 

We shortly discuss presently running and upcoming experiments. Basically all 
of them will use enriched material, and all are located in underground laboratories. 
GERDA34 ^jj^ MAJORANA^^ wiU use ^'^Ge, in the case of GERDA operated in 
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liquid Argon. Phase I consists of 18 kg previously used by IGEX and Heidelberg- 
Moscow and will test in the very near future (till 2013) the Heidelberg-Moscow 
claim, which unambiguously will only be possible with the same isotope. Phase II 
will work with 40 kg and depending on the outcome, a phase HI, probably joined 
with MAJORANA (2 x 60 kg), is possible. Inauguration of GERDA was in Novem- 
ber 2010. CUORE^^ extends CUORICINO to several towers of material, aiming at 
200 kg "OTe and a start of data taking in 2013. EXO^^, whose prototype with 200 
kg of liquid Xenon enriched to 80% is in commissioning, will apply liquid or gaseous 
Xenon; by using a time-projection chamber there is sensitivity to event topology. It 
will be attempted in a later phase to laser-tag the ^'^^Ba"'""'" ion, which is the decay 
product of the isotope under investigation, -^'^^Xe. SuperNEMO^-^ uses the NEMO-3 
approach and wiU work with about 100 kg of ^^Se or ^^°Nd. SNO-|-^° wishes to fill 
the large SNO detector with a total of 44 kg of ^^^Nd. KamLAND-Zen^'^ pursues 
a similar approach with the KamLAND experiment, using ^'^^Xe. CANDLES^^ 
will investigate CaF2 scintillators, and is currently analyzing enrichment options 
for later phases. COBRA'^° will be an array of CdZnTe room temperature semi- 
conductors, mainly sensitive to Oi^/3/3 of ^^^Cd, but to other decay modes as well. 
LUCIFER'^^ proposes to use scintillating bolometers at low temperature. MOON'^® 
wants to use scintillators in between source foils, DCBA'^^ aims at putting source 
foils with i5°Nd in a magnetized drift chamber. XMASS proposes liquid scintil- 
lating Xenon, NEXT'^^ a gaseous Xenon TPC. Some of the experiments can also 
be used as solar neutrino or dark matter experiments, such as XMASS, NEXT or 
MOON. More details on the experiments can be found in the respective publications 
and the reviews®'^'^^'^^. 

It is encouraging that different experimental techniques will be pursued, and 
that different isotopes are under study. Eventually, a multi-isotope determination 
of Oiy/Sp would be preferable, to make it more unlikely that a peak coming from 
an unidentified background process mimics the signal. This is the first reason for 
multi-isotope determination. 

We will focus in this review on neutrino-less double beta decay. However, there 
are similar processes called neutrino- less double beta"*" decay (Oj^/3"'"/3"'"), or beta"''- 
decay electron capture (Oi//?+EC), or double electron capture (Oi^ECEC) of bound 
state electrons e^, which can also be searched for: 



Observation of one of those processes would also imply the non-conservation of lep- 
ton number. The rates depend on the particle physics parameters in the same way 
as Oi//3(3 does. The creation of two positrons reduces the phase space and renders 
rates for 0i//3~^/3~^ very low. Somewhat less suppressed are (Oi^/3"'"EC) processes. In 



{A,Z) ^ {A,Z -2) + 2e+ 
+ {A,Z)^ {A,Z-2) + e+ 
2e^ + {A,Z) ^ {A,Z-2y 



(Ot^/3+EC) , 
(Oi/ECEC) . 



(7) 
(8) 
(9) 
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Oi'ECEC the final atom (and sometimes the nucleus) are in excited states and gen- 
erate photons (and 7 rays). The rate is low, unless a resonance can be met. This 
occurs^^ for certain Oz^ECEC modes, if the initial and final states of the system 
are degenerate in energy. Here the ^^^Gd-^^^Sm transition has, via Penning-trap 
mass-ratio measurements, recently been identified as an interesting candidate for 
neutrino- less double electron capture^^, though it is currently unclear if an ex- 
periment competitive to Oj^/3/3-searches can be realized. The current limits of the 
reactions (7,8,9) are summarized in^^. The main focus of future experiments is on 
the standard process Qv^jS in Eq. (1). The other reactions could however be used 
to distinguish different Oi^/J/J-mechanisms from each other, see Section 6. 

3. Nuclear physics aspects 

Nuclear physics is (unfortunately) an almost irreducible difficulty in making inter- 
pretations of neutrino-less double beta decay. Observation of the process means 
of course the proof of lepton number violation, but more precise particle physics 
interpretations suffer from any nuclear physics uncertainty. The calculation of the 
Nuclear Matrix Element (NME) is a complicated many body nuclear physics 
problem as old as Qvfiji. It basically describes the overlap of the nuclear wave 
functions of the initial and final states. A nuclear model typically has a set of 
single-particle states with a number of possible wave function configurations, and 
diagonalizes a Hamiltonian in a mean background field. A general property of solv- 
ing Hamiltonians is that the energy levels are rather stable in what regards small 
modifications. Wave functions, and hence overlap, are however very sensitive to 
small modifications of the Hamiltonian, and this is the origin of the uncertainty in 
the values of NMEs. 

3.1. Standard mechanism 

Most theoretical work has been invested into the study of the standard mechanism 
of light neutrino exchange, on which we will focus in the following discussion. The 
process is evaluated as two pointlike Fermi vertices and the exchange of a light 
neutrino with momentum of about q ~ 0.1 GeV, corresponding to the average 
distance r~ l/? — 1 fm between the two decaying nuclei. Since the neutrino is 
very light with respect to the energy scale one denotes the situation as a "long-range 
process" . 

The expression for the decay rate is 

F""^ = G"''(Q,Z)|Al"n'-^^^, (10) 

mi 

with the phase space factor G°'^((5,2'), (mee) the particle physics parameter'^ in 
case of light neutrino exchange (to be defined in Section 4), and M.^'^ the NME. 



'^Sometimes one includes the electron mass in the phase space factor. 
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The quantity {niee) is usually called the "cfFcctivc mass", or the "effective elec- 
tron neutrino mass" . The 5 main approaches to tackle the problem arc the Quasi- 
particle Random Phase Approximation (QRPA, including its many variants and 
evolution steps)*^'*^, the Nuclear Shell Model (NSM)^^, the Interacting Boson 
Model (IBM)'™, the Generating Coordinate Method (GCM)^^ and the projected 
Hartree-Fock-Bogoliubov model (pHFB)^^. We will not go into comparing in de- 
tail the different procedures, and refer the reader to the cited papers and the 
reviews^''^'^'^'^'^'®'^'^^. One example on how the approaches differ is to note that 
QRPA calculations can take into account a huge number of single particle states 
but only a limited set of configurations, whereas in the NSM the situation is essen- 
tially the opposite. The issue of which method should be used is far from settled. 
Another point arc short-range correlations (SRC), since the main contribution to 
NMEs comes from internucleon distances r ^ (2 — 3) fm^'^, and the nucleons tend to 
overlap. SRC take the hard core repulsion into account. There are different proposals 
on how to treat SRC, namely via a Jastrow-like function. Unitary Correlation Op- 
erator Method (UCOM), or Coupled Cluster Method (CCM). The Jastrow method 
leads typically to a reduction of NMEs by about 20% while UCOM and CCM both 
reduce NME by about 5% as compared to calculations without SRC^'^''*^. 

In contrast to 2j//3/3, which involves only Gamov- Teller transitions through in- 
termediate 1+ states (because of low momentum transfer), 0^/3/3 involves all mul- 
tipolarities in the intermediate odd-odd (A, Z + 1) nucleus, and contains a Fermi 
and a Gamov- Teller part (plus a negligible tensor contribution from higher order 
currents): 

The matrix elements for the final and initial states |/) and |i) can be written as 

■^GT = (/I E <^k Tf T- HGT{rik,Ea)\l) , 

= (/I E rr r,7 i/F(m, i?a)N) , ^^^^ 

Ik 

where rik — l/g — 1/(0.1 GeV) is the distance between the two decaying neutrons 
and Ea is an average energy (closure approximation due to the large momentum of 
the virtual neutrino). The "neutrino potential" 

oo 

H{x,y) oc i / dq — (13) 



integrates over the virtual neutrino momenta. The two emitted electrons are usually 
described in s-wave form because one focusses on O^g — > Og g transitions, p-wave 
emission, which would lead to transitions to excited states, is suppressed in the 
standard neutrino exchange mechanism^, see Section 6.2. The 2vj3l3 matrix elements 
can be written as (note the different energy dependence in comparison with the 
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Table 5. Dimensionless NMEs calculated in different frameworks, normalized to rg = 1.2 fm 
and QA = 1.25. The method used to take into account short range correlations is indicated in 
brackets. There is also a pseudo-5f/(3) model^* for the highly deformed nucleus ^'''"Nd, with a 
matrix element 1.00. Taken from^^. 





NSM 


Tiibingen 


Jyvaskyla 


IBM 


GCM 


pHFB 


Isotope 


(UCOM)« 


(CCM)47 


(UCOM)-** 


(Jastrow)^" 


(UCOM)5i 


(mixed)^2 


«Ca 


0.85 








2.37 




76 Ge 


2.81 


4.44 - 7.24 


4.195 - 5.355 


4.636 - 5.465 


4.6 




82 Se 


2.64 


3.85 - 6.46 


2.942 - 3.722 


3.805 - 4.412 


4.22 




96 Zr 




1.56 - 2.31 


2.764 - 3.117 




5.65 


2.24 - 3.46 


ioomo 




3.17 - 6.07 


3.103 - 3.931 


3.732 - 4.217 


5.08 


4.71 - 7.77 


llOpd 












5.33 - 8.91 


116Cd 




2.51 - 4.52 


2.996 - 3.935 




4.72 




"4Sn 


2.62 








4.81 




l30Te 


2.65 


3.19 - 5.50 


3.483 - 4.221 


3.372 - 4.059 


5.13 


2.99 - 5.12 


I36xe 


2.19 


1.71 - 3.53 


2.38 - 2.802 




4.2 




isONd 




3.45 




2.321 - 2.888 


1.71 


1.98 - 3.7 



NMEs for OujSp) 

(/I E \n) {'n\ J2 Ob \i) 

^^^^ E^^{M.-Mf)/2 ' ^ ^ 

(14) 



{f\Y.r-\n){n\Y.r^\, 

a b 

En - {Kh - Mf)/2 



where the sum over n includes only 1+ states. This is the reason why 2j//3/3 gives 
only indirect information on 

We would like to stress here that care has to be taken when different calculations 
are compared^^'^^: for instance, NMEs are made dimensionless by putting a factor 
1/R\ = l/{roA3)'^ in the phase space factor (in the convention of Eq. (11) the 
phase space becomes independent of g^i), where in the nuclear radius Ra the pa- 
rameter tq is sometimes chosen as 1.1 fm or 1.2 fm. The axial- vector coupling gA is 
often chosen to be 1.25 or 1.0. In addition, it is often overlooked (see the discussion 
jj^54,56-j ^Yia.t the phase space factors G''^'^{Q,Z) can differ by up to order 10%, for 
instance when one compares the results from^ or^^. The results from^ are given in 
Table 1. In Table 5 and Fig. 3 we give a compilation^^ of NME values from different 
calculations (see Refs.^^'®°'^^ for similar recent compilations). For definiteness, we 
will often apply the values of this table in what follows. Main features of the current 
status are that NMEs using QRPA seem to agree with each other and also with 
IBM calculations. NSM evaluations are consistently smaller and show little depen- 
dence on Z . However, it is encouraging that conceptually different approaches give 
results in the same ballpark. All in all, there has been some improvement in recent 
years, in particular the number of groups and approaches, as well as experimental 
support, has increased. However, full understanding and/or consensus has not been 
reached yet. Currently, one has to take the uncertainty at face value and keep it in 
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QRPA (Jy) - ■ - ■ 

1 IBM 

T ■ IBM 
T ; I GCM ♦ 

I I PHFB - ■ - 

g _ I T ^ I . Pseuclo-Sy(3j V 

> L I * I I : T 

» ± \± - i- ^ . I » 



-■'ft: h t 

T . 11 _ I 

T a 

I 
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Fig. 3. Nuclear Matrix Elements for Oi^/3/3, different isotopes and calculational approaches. 'Tue' 
and 'Jy' are both QRPA results. 



mind when interpreting the results of OvjSp experiments. Table 2 gives the current 
limits on the effective mass (mge) obtained with the NME compilation from Table 
5 and the phase space factors from Table 1. The best limit on {rriee) is provided 
by ''^Gc. but, as we will see, for some other mechanisms stronger limits stem from 
other isotopes, in particular ^■^°Te. 

If available, we include uncertainties of the NMEs in Table 5 and Fig. 3. How- 
ever, not all authors provide errors in their calculations. Those theoretical uncer- 
tainties can arise from varying gA, gpp (the particle-particle strength parameter in 
QRPA models), or other model details. One can distinguish here between correlated 
(e.g. the use of SRC or the value of gA or gpp) and uncorrelated errors (e.g. the 
model space of the single particle base)^^. Eventually, a multi- isotope determina- 
tion of 0iyf3/3 would be preferable, to disentangle the different types of errors. This 
is the second reason for multi-isotope determination. Ideally, if there was 
one adjustable parameter x in the calculations, then two measurements would suf- 
fice to fix (rriee) and x (up to degeneracies). A third result would overconstrain 
the system^'^ and allow for cross checks. This requires analyses of degeneracies and 
realistic estimates of theoretical errors, an effort which has recently started^^. With 
the factorization in Eq. (4) it is clear that the ratios of two measured half-lifes are 

T°J,{Ai,Zi) ^ GiQ2,Z2)\MiA2,Z2)\^ 

Tl-^{A2,Z2) G{Qi,Zi)\M{Ai,Zi)\^' ^ ' 

i.e. the particle physics parameter drops out. The ratio is sensitive to the NME 
calculation^^, and systematic errors are expected to cancel. Fig. 4, taken from^^, 
shows the error ellipses of matrix elements within a QRPA analysis. 
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Free parameters of an Ansatz can in principle also be fixed or tested by other 
means, for instance the particle-particle strength parameter gpp in QRPA models 
can be adjusted to reproduce the 21/(5(3^^, single beta decay^^ or electron capture 
rates. Overconstraining the parameters is possible if data on all these processes is 
available^^. 

In recent years the uncertainty of the individual approaches to NMEs, which 
was somewhat overestimated in the past^^, has been reduced. An experimental 
program to support the calculations with as much information as possible was 
launched^^, including charge exchange reactions^^ to determine Gamov- Teller tran- 
sition strengths. The latter are directly related to 2i/^/3 matrix elements; applying 
the results to Oi//?/3 requires theoretical input, see e.g."^". Occupation numbers of 
neutron valence orbits in the initial and final nuclei are not known very well, and 
measurements^"'^ via nucleon transfer reactions are helpful for all NME approaches. 
Muon capture rates can also be usefuF^ , because the momentum transfer is of order 

~ 100 MeV, i.e. of the same order as for Qv^p. Determinations of Q-values^^ 
with precision spectroscopy^^ are also ongoing. This is particularly helpful for ex- 
periments in which the energy resolution is comparable to the current uncertainty 
of the Q-value. Another motivation for precise Q- value determinations is the iden- 
tification of candidates for Oz^ECEC which show resonance behavior, as mentioned 
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Fig. 4. 1(7 error ellipses (logarithms of the NMEs) within QRPA calculations. The major axis 
corresponds to variation of the short-range correlation model (blue is Jastrow, red is UCOM) and 
qa, while the minor axis corresponds to variations of gpp. Taken from^^. 
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above. 

The QRPA particle-particle strength parameter gpp can be fixed by the mea- 
sured 2^/3/3 rates and used as input for Oi^/3/3 predictions. However, using this gpp 
value for calculating the rates of beta-decay or electron capture of the intermediate 
double beta decay isotope sometimes fails. One particular observational approach 
to these issues is the TITAN-EC experiment''^, which aims at testing with ion traps 
the badly known electron capture rates of the intermediate odd-odd state of double 
beta decay via observing the de-excitation X-rays. 

It is expected that the future the uncertainty in the NMEs will further decrease, 
though (owing to the enormous complexity of the problem) the precision by the end 
of the decade will probably not be better than 20%. 

3.2. Non-standard mechanisms 

The evaluation of NMEs in non-standard mechanisms is a less well developed field, 
with less calculations available, and often only within one particular nuclear physics 
approach. In general the NMEs can obtain now contributions from Fermi, Gamov- 
Teller, pseudo-scalar, tensor, etc. contributions, and the realization of Oj^/3/3 can 
differ from the standard mechanism in 

(i) the Lorentz structure of the currents (e.g. right-handed currents); 

(ii) the mass scale of the exchanged particle (e.g. exchange of heavy SUSY parti- 
cles); 

(iii) the number of particles in the final states (e.g. modes with additional Majoron 
emission); 

A frequent feature here is that the scale of Icpton number violating physics is 
larger than the momentum transfer or nuclear energies, in which case one speaks of 
a "short-range process" . Non-standard physics including light neutrino exchange, 
hence long-range, is however also possible. Within QRPA, a general Lorentz- 
invariant parametrization of the Qv^fi decay rate has been developed for long-^^ 
and short-range^^ processes (these papers are in fact the only entries on SPIRES 
with the word "superformula" in the title). In those works the most general La- 
grangian for Oj^/3/3 was written down and each term includes an individual prefactor 
Ci. These can in principle via Fierz-transformations be translated (see also''®) into 
the particle physics parameters of the alternative realizations of Oj^/3/3 which we will 
discuss in Section 5. 

It is possible that the different Lorentz structure leads to additional contribu- 
tions to the NMEs, which are not present in other realizations. Different Lorentz 
structure implies also that the energies of the individual electrons, and their angu- 
lar distribution, may differ from the standard mechanism'^' ''^'®°. Further potential 
differences are modified relations between the rate of 0^/3/3 and Oi'/J+EC®^, or with 
the decay rate to excited states®^'®'^'®^'®^. Some details will be discussed in Section 
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6.2. Short distance physics implies that the heavy particles with mass Mx can be 
integrated out and the pure particle physics amplitude is inversely proportional to 
Mx or M^, depending on whether it is a fermion or boson. A nuclear aspect of 
short distance physics is that the inner structure of the nucleons becomes relevant, 
which is taken into account by multiplying the weak nucleon vertices with (dipole) 
form factors^^ 

9A 

(1 - q'/MlY 

with a mass parameter M\ ~ (0.9 GeV)^. This introduces e.g. for heavy neutrino 
exchange a dependence proportional to M\, after an integration over the neutrino 
momenta has been performed in the potential Eq. (13). Here the form factor avoids 
the otherwise exponential suppression of the amplitude due to the repulsion of the 
nuclei. Finally, if additional particles are emitted in addition to the two electrons, 
such as in Majoron modes (Section 5.5), significant phase space effects can be 
expected. 

Another aspect of heavy particle exchange is that pion exchange can 
dominate^^'®^. This means that the pions which are present in the nuclear medium 
undergo transitions like tt^ — > 7r+e^e~, i.e. the hadronization procedure of the 
quark level diagram differs from the 2 nucleon mode discussed so far. Though the 
probability of finding pions in the nuclear soup is less than 1, this can be compen- 
sated by the fact that the suppression due to the short-range nature of the usual 2 
nucleon mode is absent, because low mass pions can mediate between more distant 
nucleons. One or both of the two initial quarks can be placed into a pion. In fact, 
i?-parity violating SUSY contributions (Section 5.4), turn out to be dominated by 
pion NMEs. 

We will discuss aspects relevant to particular non-standard mechanisms and 
means to distinguish them, from one another and from the standard one, in the 
later sections which deal with the respective mechanisms. Different realizations of 
the decay influence the NMEs in a way which depends on the isotope and on parti- 
cle physics. Therefore, eventually a multi-isotope determination of Oz^/3/3 would be 
preferable, in order to disentangle the different mechanisms^^'^"'^^'^^'^"^, see Section 
6. This is the third reason for multi-isotope determination. 



4. Standard Interpretation 

In this Section we will discuss the standard mechanism of neutrino-less double beta 
decay, let us repeat for convenience the definition: 

Neutrino-less double beta decay is mediated by light and massive Majorana neu- 
trinos (the ones which oscillate) and all other mechanisms potentially leading to 
0yj3j3 give negligible or no contribution. 

We will first summarize the current status of our understanding of lepton mix- 
ing and neutrino mass, before discussing the amount of information encoded in 
01^/3/3 combined with the standard interpretation. Readers who are very familiar 
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with neutrino physics can go directly to Section 4.2 and skip the summary of neu- 
trino physics in Section 4.1. 

4.1. Neutrino physics 

Most part of the review will deal with the standard, and presumably best motivated, 
interpretation of QvjSfi, light massive Majorana neutrino exchange. We will first 
review the current theoretical and phenomenological status of neutrino physics. 

4.1.1. Neutrino mass and mixing: theoretical origin 

The theory behind neutrino mass and lepton mixing has been reviewed in several 
places^^'^^'^^'^^. Lepton mixing is rather different from quark mixing. In addition, 
the mass of the two lepton partners in an SU (2) ^ doublet (e.g. i>e and e) is extremely 
hierarchical, in sharp contrast to the partners in quark doublets (e.g. u and d, 
with ruu = 0(md))- It is natural to believe that these discrepancies are related to 
special properties of the neutrinos. Indeed, most, if not all, theorists believe that 
neutrinos are Majorana particles. This is the case in basically all Grand Unified 
Theories (GUTs), and also from an effective theory point of view, in which non- 
renormalizable higher dimensional operators invariant under the SM gauge group 
are constructed. The lowest dimensional (Weinberg) operator is unique, and reads^^ 

2 A " — (17) 

Here the superscript denotes the charge-conjugated spinor. La = {vonOcY' are 
the lepton doublets of flavor a = e,/x,T and $ is the Higgs doublet with vacuum 
expectation value w = 174 GeV. A Majorana neutrino mass matrix is induced by 
this operator, given by m^, = hv^ /h.. With the typical mass scale of ~ 0.05 eV, 
it follows that A ~ 10^^ GeV, tantalizingly close to the GUT scale. This is one of 
the main reasons why neutrino physics is popular: large scales are probed by small 
neutrino masses. It has been shown that within the minimal standard electroweak 
gauge model, there are only three tree-level realizations®^ of the Weinberg operator. 
One is the canonical type I seesaw mechanism with right-handed neutrinos. An- 
other approach is introducing a scalar Higgs triplet (type II, or triplet seesaw^*^^), 
and the third one involves hypercharge-less fermion triplets (type HI seesaw^^^). In 
Table 6 wc summarize the main approaches for generating small neutrino mass. 

Taking first the standard type I seesaw as an example, one introduces 3 (actually 
2 would suffice) Majorana neutrinos A^ii.i, which have a Majorana mass matrix M^. 
After electroweak symmetry breaking a Dirac mass term with the SM neutrinos is 
present and the full Lagrangian for neutrino masses is 

'^ = ~\M-^^^^'.-^oir,., = -\T^^ (18) 

with Nn = {Nri, Nr2, Njis) and i^l ~ (i^e, i^/^, ^'t)l- We see that the combination 
of a Dirac and a Majorana mass term is a Majorana mass term, no matter how 
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Table 6. Tree-level approaches to small neutrino mass classified according to the ingredient which has to be added to the SM, 
and the electroweak quantum numbers of the new particles. LNV denotes Lepton Number Violation. 
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hSL^ + TrSMsS 


(hvr 

Ms 


A = 1 A/j: 



small the Majorana mass is. Being SM singlets, the scale of Mr is not connected to 
the only energy scale of the SM (the Higgs vacuum expectation value), and hence 
can be arbitrarily high. Integrating out the heavy states, or block-diagonalizing'^ 
the mass matrix in Eq. (18) gives a Majorana mass term for the light neutrinos, 

= ~mj) Mj^ ^ niD , (19) 

plus terms of order m^/M^. The states for which this mass matrix is valid are the 
initial vl plus a contribution of N'^, which is however suppressed by mo /Mr. We 
see that the Weinberg operator is realized with A ~ Mr. 

Often one considers a triplet term for neutrino masses, generated by an SU(2) l 
triplet scalar with non-zero vev vl of its neutral component. The coupling of the 
triplet to two lepton doublets with the Yukawa coupling matrix h gives a neutrino 
mass rriu = Ml, i.e. a direct contribution (see Section 5.2). Of course, both the 
type I and the type II term could be present. In this case the zero in the upper left 
entry of Eq. (18) is filled with a term M^. The neutrino mass matrix in this case 
reads 

m,, = Ml- Tn^M~^ mo ■ (20) 

Finally, type III seesaw introduces 3 hypcrchargc-lcss fermion triplets (one for each 
massive light neutrino), whose neutral components play the role of the Nm of type 
I seesaw. 

What about production of low scale seesaw messengers at colliders? A recent 
review on the situation can be found in^*^"^. While Majorana neutrino production 
proves difficult because of the constraint of its small mixing with SM particles, 
Higgs and fermion triplets have gauge quantum numbers and can be observed up 
to TeV masses. Note that in left-right symmetric models, or models with gauged 



''The condition for this is that the eigenvalues of AIr arc much heavier than the entries of m^j. 



January 27, 2013 17:31 WSPC/INSTRUCTION FILE DoubleBetaReview 



20 W. Rodejohann 

B — L, the right-handed neutrinos can have gauge interactions and can be produced 
more easily at colhders. The lepton number violation associated with the seesaw 
messengers can lead to their identification and spectacular like-sign lepton events. 
For this to be realized one needs to bring the seesaw scale down to TeV, see^°* for 
a recent review on how to achieve this. 

It is clear that, either way, neutrinos are Majorana particles, i.e. 

„^ = CdT = u,. (21) 

Here we have chosen a convention in which there is no phase in the above relation. 

In general the mass matrices for neutrinos (m^) and for charged leptons (me) 
are non-trivial. Diagonalizing those matrices with unitary'^ Ui, and Ui, respectively, 
results in the charged current term in the appearance of the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix U = Uj [/^: 

Ccc = -^ 4 7^ Uc.^ V, W- . (22) 
v2 

In the basis in which the charged leptons are real and diagonal, the neutrino mass 
matrix is diagonalized by U . It is useful for our purposes to stay in this basis, in 
which the mass matrix for Majorana neutrinos can be written as 

= U* TO^J'^'s fjt ^ ^here to^^'^s = diag(mi, ms, mg) . (23) 

The mass matrix is complex and symmetric; after rephasing of three phases there are 
9 physical parameters. Because the mass term goes as ^^'^ cx v, the Lagrangian 
is not invariant under a global transformation f — )• e^'^i'. The charge associated 
with this transformation, lepton number, is therefore not a conserved quantity and 
L is violated by two units'^. This is exactly what is required for the presence of 
neutrino-less double beta decay. 

Another appealing prediction of seesaw is the possible generation of the baryon 
asymmetry of the Universe via leptogenesis^"^. Here the heavy seesaw messengers 
decay out of equilibrium (Sakharov condition I) in the early Universe and, due 
to CP violating phases (condition II), create a lepton asymmetry which subse- 
quently is transfered into a baryon asymmetry via B + L violating (condition HI) 
non-perturbative SM processes. In this context, proving the Majorana nature of 
neutrinos and the presence of CP violation in the lepton sector would strengthen 
our belief in this already very appealing mechanism. This remains true even though 
a model-independent connection between the low energy CP phases and the neces- 
sary CP violation for leptogenesis cannot be established^"^. In general, taking the 

"^Strictly speaking the matrix Ui, is not unitary in type I seesaw, due to mixing of the leptons 
with the heavy neutrinos. This is however usually a very small effect \U,^Ul — 1| ~ {mo/Mn)^ 
and phenomenologieally constrained to be less than a permille effect^''^ . 

"^It should be noted that there are alternatives to the seesaw mechanism, see^"® for a discussion. 
Examples are radiative mechanisms, supersymmetric scenarios, or extra-dimensional approaches. 
It can happen that lepton number is conserved in such frameworks and Ouf}^ cannot take place, 
but this is clearly a rare exception rather than the rule. 
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standard type I seesaw as an example, there are in total six CP phases, three of 
which get lost when the heavy Majorana neutrino mass matrix is integrated out 
to obtain rrii, (see Eq. (19)). In principle, one could construct models in which the 
"low energy phases" take CP conserving values, while the remaining three phases 
are responsible for leptogenesis®. Leaving this seemingly unnatural possibility aside, 
one expects that CP violation in the lepton sector at low energy is present if there 
is "high energy" CP violation responsible for leptogenesis. One should however not 
expect that the Majorana phases are "more connected" to leptogenesis than the 
Dirac phase. At the fundamental (seesaw) scale, there are six CP phases and the 
three low energy phases will be some complicated function of these phases and 
the other seesaw parameters. From this point of view, the low energy Dirac and 
Majorana phases are not different from each other. 

A final remark necessary here is that a link between QvjSp and the baryon 
asymmetry is not guaranteed. The often-made and popular statement that Qv^jS- 
experiments probe the origin of matter in the Universe is not true. For instance, if 
neutrino mass is simply generated by a Higgs triplet, then this triplet alone cannot 
generate a baryon asymmetry, but QvjSp is very well possible. 

4.1.2. Neutrino mass and mixing: observational status 

Neutrino oscillations have been observed with solar, atmospheric and man-made 
(reactor, accelerator) neutrinos, see^"^'^*^^ for extensive reviews on the status of 
neutrino physics. This implies that in the charged current term of electroweak 
interactions the neutrino flavor states Ve, and Vr are superpositions of neutrino 
mass states: 



where a — e,/i, r and i = 1,2,3. The PMNS mixing matrix U is unitary and can 
be written in its standard parametrization as 



where s^- = sin 9ij , Cij = cos 9ij and 5 is the "Dirac phase" responsible for CP viola- 
tion in neutrino oscillation experiments. This phase is expressible in a parametriza- 
tion independent form as a Jarlskog invariant: 

Jcp = Im {C/*i C/e3 C^Mi} = ^ sin 2^12 sin 26*23 sin 26*13 cos 6*13 sinJ , (26) 

where we have given it in its explicit form for the standard parameterization. In 
Eq. (25) we have included a diagonal phase matrix P, containing the two "Majorana 

*Note that the opposite case is also possible. 
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(24) 




(25) 
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phases" a and 13: 

P = diag(l,e'",e'('^+^)). (27) 



These phases are physical^^'^ if neutrinos are Majorana particles. Note that we have 
included i5 in P, in which case the first row of the PMNS matrix is independent of 
5. For three neutrinos we have therefore 9 physical parameters, three masses toi 2,3, 
three mixing angles 6*12, 6*13, ^23 and three phases (5, a, /3. 

One can also define invariants for the Majorana phascs^^^, for instance Si = 
lm{UeiU*2} = -C12 512 0^3 sin a and S2 = Im{C/e2t^e3} = S12C13S13 sin((5 - ^). 
Note that CP violation due to the Majorana phases is present only if, in addition 
to Si = Im{J7eiC/e2} 7^ 0; Re {C/ei [7*2} ^ also holds. The reason for this is 
that the cases a, /3 = 7r/2 correspond to the CP parities of the Majorana fields^^^, 
which can be either positive or negative. Majorana phases are present because the 
mass term in the Lagrangian is proportional to {m,j)ai3V^Vi3 and a rephasing of 
the spinors Va can eliminate fewer phases than in the Dirac case, where the mass 
term is (TOi/)a/3 Va vp- For N Majorana neutrinos, there are N — 1 Majorana phases 
in addition to ^ {N — 2) {N — 1) Dirac phases and [N — 1) mixing angles^. 

For three fermion families, neutrino oscillation experiments are sensitive to the 
three mixing angles, the two independent mass-squared differences (including their 



^This counting is valid for active neutrinos only. For N massive families including 7^ Ns = N — 3 
massive sterile neutrinos, one has A'^ — 1 = Ns + 2 Majorana phases, 3 (A'^ — 2) = 3 {Ns + 1) mixing 
angles and 2 TV — 5 = 2 Ns + 1 Dirac phases. The number of angles and Dirac phases is less because 
the i TVs (TVs — 1) rotations between sterile states are unphysical. 
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Table 7. Current values from global fits to the world's neutrino os- 



dilation experiments, 
the inverted ordering. 


Taken from^^^ 


The values in 


brackets are for 


parameter 


best-fit j^j^^ 


2a 


3(7 


Am|i [10-5 eV^] 


7 fi4+0-i9 
'■"^-0.18 


7.27 - 8.03 


7.12 - 8.23 


|Am2j [10-3 eV^] 




2.28 - 2.64 
(2.17 - 2.54) 


2.18 - 2.73 
(2.08 - 2.64) 


sin^ 012 


0-316Io.016 


0.29 - 0.35 


0.27 - 0.37 


sin^ 6*23 


O.Sliro^e 
(0.52lH«6) 


0.41 - 0.61 
(0.42 - 0.61) 


0.39 - 0.64 
(0.39 - 0.64) 


sin^ 6*13 


o.oi7±n;!^ 


< 0.031 
(< 0.036) 


< 0.040 
(< 0.044) 



sign), and the Dirac phase 5. The general formula for oscillation probabilities is 
P{i^c. ^ ^^/j) = - 4 E Re {C/*, U^^ Up, C/„, } sin^ 

>, A 2 r 

C/*,C/;^.C/^,C/„4sin^^, 

with E the neutrino energy and L the baseline. To leading order, using the hierarchy 
of AttIq = Ato|]^ <C IAtoj]^! ~ I AtoIjI = Am^, this formula usually breaks down to 
two neutrino oscillation formulas. The angle 9i2 and ^ni^i = Amg are responsible 
for solar neutrino (suitably modified with matter effects) and long-bascline reactor 
neutrino oscillations. Atmospheric neutrinos are governed by 023 and Am^j, the 
same parameters which long-baseline accelerator neutrinos are sensitive to. Finally, 
013 (if non-zero) and Am|i are responsible for short-baseline reactor neutrino and 
long-baseline — > J^e oscillations. Non-zero also provides a link between the so- 
lar and atmospheric sector and is intensively searched for, as leptonic CP violation 
in oscillations would be absent if it was zero. Our current knowledge of the oscilla- 
tion quantities is summarized in Fig. 5 and Table 7, taken from^^'^. It is noteworthy 
that the sign of the (atmospheric) mass-squared difference is unknown, as are the 
three CP phases, thus including the Majorana phases. The hint^^^ towards non- 
zero 013 recently exceeded the 3a level^^^, after the T2K long-baseline experiment 
provided evidence for electron neutrino appearance^^^. 

A recent review on the details of current and future determinations of the pa- 
rameters can be found in Ref.^*^®. An extensive program to improve the precision on 
^13, ^23 and Am|i (including its sign) has been launched, while improvement in the 
precision of 0i2 and Am^i is not on top of the neutrino community's agenda, mostly 
because they are currently the best-known parameters. As we will see below, pre- 
cision determination of the solar parameters may be required if the inverted mass 
ordering is to be tested with Qv[ij3. The unknown sign of the atmospheric mass- 
squared difference defines the mass ordering (see Fig. 6): normal for Am^ > 0, 
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=3 
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■ V, i 







V2 
Vi 



V3 



Fig. 6. Normal (loft) vs. inverted (right) mass ordering. The red area denotes the electron content 
|J7eiP in the mass state Vi. Accordingly, the yellow and blue areas denote the muon and tau 
contents. Taken from-*^"^. 



inverted for Am^ < 0. The two larger masses for each ordering are given in terms 
of the smallest mass and the mass squared differences as 

normal: TO2 
inverted 



Ato| 



TO2 



Ato| 



Ami ; 



nil ~ yml + Ar 



(29) 



Note that the oscillation data and the possible mass spectra and orderings are 
independent on whether neutrinos are Dirac or Majorana particles. The two possible 
mass orderings arc shown in Fig. 6. Of special interest are the following three 
extreme cases: 



normal hierarchy (NH): 7713 
inverted hierarchy (IH): 
quasi-dcgencracy (QD): 



a/AtoJ > 7712 — \ AniQ > 7771 



7772 

2 



7771 

-i2 ^ 



x/Ami > 7773 , 

.2 ^ ^2 ^ A^2 



(30) 



Q = nil — 11L2 — 7773 ^ AtTI^ , 

As can be seen from Fig. 5 and Table 7, the current data is well described by so- 
called tri-bimaximal mixing^^'', corresponding to sin^ 6I13 = X cos^ 6'i3, sin 612 = 
i X cos^ 612 and sin^ 6*23 = 1 x cos^ 623 : 



U 







V 



(31) 



The application of flavor symmetries to the fermion sector, in order to obtain this 
and other possible mixing schemes is a very active field of research. For references 
and an overview of flavor symmetry models, see^^^. 



A longstanding issue in oscillation physics is the indication of the presence of 
sterile neutrinos. The LSND expcrimcnt^^^ found evidence for 77^ — )■ Ve transitions 
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which, when interpreted in terms of oscillations, are described by a Am^ ^ eV^ and 
small mixing < 0.1. These values survive even when combined"'^^^'"'^^^'^^^ with the 
negative results from the KARMEN experiment^^^. This mass scale cannot be com- 
patible with solar and atmospheric oscillation and hence a fourth, sterile neutrino 
needs to be introduced. So-called 1-1-3 (3-1-1) scenarios would then be realized, in 
which one sterile neutrino is heavier (lighter) than the three active ones, separated 
by a mass gap of order eV. The MiniBooNE experiment was designed to test the 
LSND scale with different L and E, but very similar L/E. The results^^'' could 
not rule out the LSND parameters, and are also compatible with the presence of 2 
sterile neutrinos^^^. In fact, the difference between MiniBooNE's neutrino and anti- 
neutrino results can be explained by two additional eV-like Am^ plus CP violation. 
Here one could envisage 2-1-3 or 3-1-2 scenarios, in which 2 sterile neutrinos lie above 
or below the three active ones, or 1+3-1-1 scenarios^^^, in which one sterile neutrino 
is heavier than the three active ones and the other sterile neutrino is lighter. Re- 
cently, reactor neutrino fluxes have been re-evaluated and an underestimation of 
3% with respect to previous results has been found^^''. The null results of previous 
very short-baseline reactor experiments can now be interpreted as in fact being a 
deficit of neutrinos, which again is compatible with oscillations corresponding to 
Am^ ^ eV^ and small mixing ^ 0.1. A recent analysis of short-baseline neutrino 
oscillation data in a framework with one or two sterile neutrinos can be found in 
Ref.^^^. The global fit improves considerably when the existence of two sterile neu- 
trinos is assumed. 

We will discuss the situation on neutrino mass from now on. Neutrino mass can 
be measured in three and complementary different ways^: 

1) Kurie-plot experiments, 

in which the non-zero neutrino mass influences the energy distribution of electrons 
in beta decays close to the kinematical endpoint of the spectrum. As long as the 
energy resolution is larger than the mass splitting, the spectrum is described by a 
function 

{E, - Q) ^{E, - Qf - mj , (32) 
where the observable neutrino mass parameter is 

m^^y^^|[/e,Pm|- (33) 

The current limit to this quantity from spectrometer approaches is 2.3 eV at 95% 
C.L., obtained from the Mainz^^o ^^^^ Troitsk"^ collaborations. The KATRIN 
experiment^^^'^'^^ has a design sensitivity of = 0.2 eV at (90% C.L.) and a 
discovery potential of mp = 0.35 eV with 5ct significance. It represents the ultimate 

^Alternatives such as time-of-fiight measurements of supernova neutrinos cannot give comparable 
limits. Other ideas^'^^ are presumably not realizable. 
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spectrometer experiment for neutrino mass, in which an external source of beta 
emitters (tritium) is used. Further improvement of the hmits must e.g. come from 
calorimeter approaches, where the source is identical to the detector. The MARE^'^^ 
proposal will use ^® ''Re modular crystal bolometers. A history of neutrino mass lim- 
its from beta decays and reviews of upcoming experiments can be found in^^^. A 
different Ansatz called Project 8 aims to detect the coherent cyclotron radiation 
emitted by mildly relativistic electrons (like those in tritium decay) in a magnetic 
field. The relativistic shift of the cyclotron frequency allows to extract the electron 
energy from the emitted radiation^'^^. In principle, MARE and Project 8 can reach 
limits of 0.1 eV. Investigation of beta spectra is usually considered to be the least 
model-dependent Ansatz to probe neutrino mass. For instance, Refs.^^'' have shown 
that admixture of right-handed currents can be not more than a 10% effect in KA- 
TRIN's determination of mp; 

2) Cosmological and astrophysical observations 

are sensitive to neutrino mass, see^'^® for a review. In particular, effects of neutrinos 
in cosmic structure formation are used to extract limits on neutrino masses. The 
quantity which is constrained by such efforts is 



familiar from the contribution of neutrinos to hot dark matter'', ^l^, = 
E/(94.57eV). Finite neutrino masses suppress the matter power spectrum on scales 
smaller than the free-streaming scale fcps — 0.8 hnii/cV Mpc~^. However, neutrino 
mass is highly degenerate with other cosmological parameters, for instance^^*' with 
the dark energy equation of state parameter uj, so that one needs to break the de- 
generacies with different and complementary data sets. Besides cosmic microwave 
background (CMB) experiments, one can use the Hubble constant (HO) measure- 
ments, high-redshift Type-I supernovae (SN) results, information from large scale 
structure (LSS) surveys, the LSS matter power spectrum (LSSPS) and baryon 
acoustic oscillations (BAO). The impact on the neutrino mass limit is shown in 
Table 8, taken from'^*'^, in which a fit to a cosmological model allowing for neu- 
trino mass, non- vanishing curvature, dark energy with equation of state uj =/= —1, 
and the presence of new particle physics whose effect on the present cosmological 
observations can be parameterized in terms of additional relativistic degrees of free- 
dom AA^ici, has been performed. As can be seen, depending on the data sets, the 
limit on S varies by a factor of 3. Future cosmological probes will add additional 
information, a summary of expectations for this is shown in Table 9. 

At the present stage it is worth noting that precision cosmology and Big Bang 
nucleosynthesis mildly favor extra radiation in the Universe beyond photons and 
ordinary neutrinos. While this could be any relativistic degree of freedom, the in- 




(34) 



''While one usually eonsiders light neutrinos as a sub-leading part of dark matter, arguments in 
favor of neutrino hot dark matter are given in^^^. 
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Tabic 8. 95% C.L. upper bound on the sum of the neutrino masses 
from different cosmological analyses. Taken from^*^. 



Model Observables rrii [eV] 



owCDM + ATVrci + 


CMB+HO+SN+BAO 


< 1.5 


ocjCDM + AA^roi + rrii. 


CMB+HO+SN+LSSPS 


< 0.76 


ACDMH 


h ruu 


CMB+HOH-SN+BAO 


< 0.61 


ACDMH 


h mv 


CMB+HO+SN+LSSPS 


< 0.36 


ACDMH 


h mv 


CMB (+SN) 


< 1-2 


ACDMH 


h my 


CMB+BAO 


< 0.75 


ACDMH 


h my 


CMB+LSSPS 


< 0.55 


ACDMH 


h my 


CMB+HO 


< 0.45 



Table 9. Future probes of neutrino mass, with their projected sensitivity. Sensitivity in 
the short term means within the next few years, while long term means by the end of the 
decade. Taken from^'^*. 



Probe 


Potential sensitivity [eV] 
(short term) 


Potential sensitivity [eV] 
(long term) 


CMB 


0.4-0.6 


0.4 


CMB with Icnsing 


0.1-0.15 


0.04 


CMB + Galaxy Distribution 


0.2 


0.05-0.1 


CMB + Lensing of Galaxies 


0.1 


0.03-0.04 


CMB + Lyman-a 


0.1-0.2 


Unknown 


CMB + Galaxy Clusters 




0.05 


CMB + 21 cm 




0.0003-0.1 




Fig. 7. 68%, 95% and 99% C.L. regions for the neutrino mass and thermally excited number 
of degrees of freedom A^s. The left plot is the A^s + 3 scheme, in which ordinary neutrinos have 
my = 0, while sterile states have a common mass scale rris, hence S ~ Nsms- The right plot is 
for the 3 + Ng scheme, where the sterile states are taken to be massless m^ = 0, and 3.046 species 
of ordinary neutrinos have a common mass my, hence S ~ 'imy. Taken from^*^. 



terpretation in terms of additional sterile neutrino species is straightforward (re- 
call the discussion on sterile neutrinos from above). Fit results very well compat- 
ible with more radiation than the SM value have been found e.g. by the WMAP 
collaboration"'^'''^ or in^'''^. This is supported by the recently reported higher ^Hc 
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abundance^"*^, which in the framework of Big Bang Nucleosynthesis can be ac- 
commodated by additional relativistic degrees of freedom, as this leads to earlier 
freeze-out of the weak reactions, resulting in a higher neutron-to-proton ratio. In 
Fig. 7 the result of a recent fit"'^^^ to cosmological data is shown, in which two situ- 
ations are analyzed: massless active neutrinos plus A^s massive sterile states {Ns + 3 
scheme); and Ng massless sterile states plus 3 massive active states (3 -f- Ng). The 
Planck satellite, with a projected sensitivity of ±0.2 to the number of extra degrees 
of freedom, will be decisive in order to test this presence of additional radiation. 
It is rather interesting that hints for the presence of sterile neutrinos are given by 
fundamentally different probes: neutrino oscillations. Big Bang Nucleosynthesis and 



Cosmological mass limits can be considered robust with respect to reasonable 
modifications of the ACDM model^^^ , in particular if different and complementary 
data sets are applied. However, several non-standard cosmologies exist for which 
no detailed study on the effect on the S bound has been performed yet, for in- 
stance coupled dark energy scenarios. Nevertheless, it is fair to say that neutrino 
masses heavier than E ~ 2 eV or so would be rather surprising, and correspond 
to very unusual scenarios. Note however that any information about the neutrino 
mass can be obtained only by way of statistical inference from the observational 
data after a parametric model has been chosen as the basis for the analysis^**®. This 
is a difference to the investigation of energy spectra in single or double beta decay 
experiments; 

3) Neutrino-less Double Beta Decay 

This possibility to test neutrino mass will be dealt with in Section 4.2.3 in some 
detail. In the ideal case, results from two or all three approaches to neutrino mass 
are present, and we will discuss this interesting case too. Neutrino mass limits 
from neutrino-less double beta decay need to assume that neutrinos are Majorana 
particles, and that no mechanism other than light neutrino exchange is responsible 
for the process. Let us note here that from Oj^/3/3 limits one can extract two different 
"masses". First, we can extract the physical masses, i.e. the eigenvalues of the 
mass matrix. These quantities are the ones tested in the other approaches. Second, 
Oi^/3/3 tests directly the quantity {mi,)ee, i-e. the ee element of the neutrino mass 
matrix in the charged lepton basis: 



Thus, the decay width is directly proportional to the fundamental quantity which 
originates at the fundamental large (seesaw) scale, without any diagonalization 
procedure dependent on known and unknown parameters. Note that in order to 
extract neutrino mass limits from 0i^/3(3 one needs to assume the neutrinos are 
Majorana particles, and that no other mechanism contributes. We will comment 
later on the complementarity of the three neutrino mass observables. 



CMB + LSS. 




(35) 
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Fig. 8. Left: quark level Feynman diagram for the standard interpretation of neutrino-less double 
beta decay. Right: geometrical visualization of the effective mass. 



4.2. Standard three neutrino picture and Ovf3f3 

In this Section wc will summarize the standard analysis of neutrino-less double 
beta decay with the standard three neutrino framework. Several works have been 
devoted in the literature to this^^^'^'**'^^^'^^'''^^^'^^2,i53,6i^ g^^^ g^rlier review can be 

found in^^^. 

The Feynman diagram for on the quark level in this interpretation is 

shown in Fig. 8. Due to the typical structure of the process it is sometimes called 
"lobster diagram" . The amplitude of the process is for the V ~ A interaction of the 
SM proportional to 

Gp 7p 1+ — 2 7- = Gp U^i— 2 7m 7+ 7<- 

^ \ " Tj2 

- Z^Gp U^, —7^ 7+ 71. , 

where 7± = ^(1 ± 75), m.^ is the neutrino mass, q ~ 100 MeV is the typical 
neutrino momentum, and Uei an element of the first row of the PMNS matrix. The 
linear dependence on the neutrino mass is expected from the requirement of a spin- 
flip, as the neutrino can be though of being emitted as a right-handed state and 
absorbed as a left-handed state. In case the interactions are not left-handed at one 
of the vertices, the linear dependence on will be absent; we will consider these 
cases later in Section 5.3. If both interactions are right-handed, the same linear 
dependence on m.i appears. Note that the amplitude is proportional to a coherent 
sum, which implies the possibility of cancellations. The decay width is proportional 
to the square of the so-called effective mass 



Im(l)| + |™e)|,2,a_^, (3), 2,/3 



(37) 



f 1 2 3^ 

which is visualized in Fig. 8 as the sum of three complex vectors mee .If one 

f 1 2 3) 

cannot form a triangle with the rriee ' , then the effective mass is non-zero. The 
Majorana phases 2a and 2/3 correspond to the relative orientation of the three 
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vectors. The standard analysis of the effective mass is the geometry of the three 
vectors expressed in terms of neutrino parameters. In the standard parametrization 
of the PMNS matrix we have 

l™ie^ I = mi pel P = mi C?2 C?3 > 

Imi'e^l -m2|C/e2p-m2S?2 4, (38) 
|m(3)| = m3|C/e3p=m3 4. 

The individual masses can, using Eq. (29), be expressed in terms of the smallest 
mass and the mass-squared differences, whose currently allowed ranges, as well as 
those of the mixing angles, are given in Table 7. From Table 2 we can read off the 
current limit on the effective mass: 



(mee)£0.5eV. (39) 

For later use we define the standard amplitude for light Majorana neutrino ex- 
change: 

A oc Gl %1 ^ 7 X 10-i« f^i) GeV-^ . (40) 
q'' \0.5 eV / 

Fig. 9 shows the future limits on the effective mass for different isotopes and half-life 
limits (see also^^'^"'^). We have again used the NME compilation from Table 5. 

The effective mass depends on 7 out of the 9 physical parameters of low en- 
ergy neutrino physics (only 623 and 6 do not appear), hence contains an enormous 
amount of information. It is the only realistic observable in which the two Majorana 
phases appear. For the other five quantities there will be complementary informa- 
tion from oscillation experiments or other experiments probing neutrino mass. It is 
also noteworthy that {rriee) is the ee element of the neutrino mass matrix m^, see 
Eq. (23), which is a fundamental object in the low energy Lagrangian. In terms of 
the origin of neutrino mass, (niee) is h^e I'^^/A, see Eq. (17) and the realizations of 
A in terms of fundamental mass scales in Table 6. 

A typical analysis of the effective mass would plot it against the smallest neu- 
trino mass, while varying the Majorana phases and/or the oscillation parameters. 
This results in Fig. 10, for which the best-fit values and 3(j ranges of the oscilla- 
tion parameters have been used. The blue shaded area is of interest because it can 
only be covered if the CP phases are non-trivial, i.e. if a,/3 7^ 0,7r/2. The values 
a,/3 = 0,7r/2 correspond to CP conserving situations, associated with positive or 
negative signs of the neutrino masses, and the resulting span of {m.ee) is also in- 
dicated in the figure. For comparison, the other mass-related observables S and 
mp are shown as a function of the smallest neutrino mass in Fig. 11. Actually, 
the smallest neutrino mass is not really an observable, so it is interesting to plot 
the effective mass against S and mp, which is shown in Fig. 12. The analytical 
expressions for the effective mass in certain extreme cases arc given in Fig. 13. 
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<n\>^^^, sin'^ 6,2 = 0.27 
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Fig. 9. Limits on the effective mass for different fialf-life limits. Tlie liorizontal lines are tlie maximal 
and minimal values of (mee) in the inverted mass ordering. 



4.2.1. Normal mass ordering 

Let us begin with the normal mass ordering. The effective mass is 



Ami si, e'^^ 



(41) 



The maximum of the effective mass is obtained when the Majorana phases are given 
hy a = (3 — 0. The case of small mi, which corresponds to a normal hierarchy (NH) 
defines the "hierarchical regime" in Fig. 13. Neglecting mi gives 



'Amlsl,cl, + ^Amlsl,e'^^--^^ , 
where both terms can be of comparable magnitude. If ^13 = one has 



mi e^*" 



Atoq + ni\ SI2 



(42) 



(43) 



where both terms can again be comparable. Note that in the last two expressions, 
as well as for other situations with small mi , the effective mass can vanish, a special 
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Fig. 10. Effective mass against the smallest neutrino mass for the 3a" ranges (top) and best-fit 
values (bottom) of the oscillation parameters. CP conserving and violating areas are indicated. 

case we will deal with in Section 4.2.6. If the smallest mass mi is much larger than 
the mass-squared differences, the effective mass for quasi-degenerate neutrinos is 
obtained: 

(mee)'^° = mo |c?2 c?3 + si, cl, e^- + sj, e'^^ \ . (44) 

Recall that mg denotes the common neutrino mass for QD neutrinos. The third 
term is now much smaller than the minimal combination of the first two terms, 
mo(ci2 — sf,), because 6*12 lies below tt/A. Therefore, the effective mass cannot 
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Fig. 12. Effective mass against sum of masses S and kinematic neutrino mass mp for the 3a" ranges 
of the oscillation parameters. CP conserving and violating areas are indicated. 

vanish for quasi-degenerate neutrinos. The estimate for the effective mass in case 
of quasi-degenerate neutrinos is 

cos 26*12 rno <, {rriee)^^ ^ mo . (45) 

This corresponds, for {rriee)'^^ ~ 0.1 eV, to half-Hfes in the regime of 10^^ to 10^^ 
yrs. Current experiments are testing this regime, thus ^ 100 kg yrs facihties with 
10^^ or less background counts are sufficient for the QD regime. 
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hierarchical 



0.1 



- 0.01 



0.001 



0.0001 



cancellation 
(only normal) 



quasi-degenerate 



0.0001 




0.001 



0.01 
m[eV] 



Fig. 13. The main properties of the effective mass as function of the smallest neutrino mass. 
We indicate the relevant formulae and the three important regimes: hierarchical, cancellation 
(only possible for normal mass ordering) and quasi-degeneracy. Cij = cos 9ij, Sij = sinOij and 
tii = tan Si, . Taken from^^'^. 



A rough estimate for the effective mass in terms of a normal hierarchy is 

, ,NH / 0.003 cV, 

(mee) ^ \ / , (46) 

W or ^Ami sin^ 9^^ < 0.003 eV) . 

The meV scale of the effective mass should be the final goal of experiments, but 
the possibility of strong or even complete cancellation has to be kept in mind. The 
half-lifes corresponding to meV effective masses are 10^® to 10^^ yrs. Multi-ton scale 
experiments are necessary for such extremely low numbers, with background levels 
below 10~^. It has been argued that if single electron events cannot be distinguished 
from double electron events, the elastic t'eC scattering of solar neutrinos represents 
an irreducible background for Oi^/3/3-experiments probing the NH regime^^^. 

In Figs. 11 and 12 we see that in case of NH lies below KATRIN's sensitivity 
of about 0.2 eV for the normal hierarchy regime. With S ~ -^/AmJ ~ 0.05 eV, 
only very optimistic or far future cosmological observations can test this value. If 
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the quasi-degenerate scenario is realized we find S ~ imp ~ 3(mee)™^'' — 3 mo- 
Table 10 shows the mass observables for the NH and QD schemes. The interplay 
of the mass observables assumes unitarity of the PMNS matrix; corrections due to 
possible non- unitarity have been discussed in Ref.^^^ and found to be negligible. 



4.2.2. Inverted mass ordering 

For the inverted mass ordering, the smallest neutrino mass is denoted ma and the 
mass matrix element is given by 



(47) 

The maximal effective mass is - as for the normal mass ordering - obtained when 
a = /? = 0. The minimal value is 



I + Ami C?2 C?3 - V mg + ^"^1 + ^'"A S?2 C?3 - "13 s?3 . (48) 



The third term of (mee) is usually negligible because ^la is small and ma is the 
smallest mass. In this case: 

(49) 



and (mee)™x = v^Am^ < (mee)™ < \/^m\ cf^ cos 26*12 = (mee) 



It is important to note that owing to the non-maximal value of 9i2 the minimal 
value of the effective mass is non-vanishing^^'^. Therefore, if limits below the minimal 
value 



= ('"-):i = (1 - iC/eaP) J Ami (1 - 2 sin' 9,,) , (50) 



' mm 



are reached by an experiment, the inverted mass ordering is ruled out if neutrinos 
are Majorana particles. If we knew by independent evidence that the mass order- 
ing is inverted (by a long-baseline oscillation experiment or a galactic supernova 
observation) then we would rule out the Majorana nature of neutrinos. Of course, 
one has to assume here that no other lepton number violating mechanism inter- 
feres. The two scales of (mee) corresponding to the minimal and maximal value of 
(mee) in case of the inverted hierarchy, given in Eq. (49), should be the intermediate 
or long-term goal of future experiments. 

The typical effective mass values of order ~ 0.03 eV are one order of magnitude 
larger than for the normal hierarchy and roughly one order of magnitude smaller 
than for quasi-degenerate neutrinos. They implies half-lifes of order 10^® to 10^^ 
yrs. A few 100 kg yrs of data taking with background levels below 10~^ counts will 
be necessary. Upcoming next generation experiments will test the inverted ordering, 
but cannot completely rule it out. 

A more detailed analysis, focussing on the important dependence on 612, will 
be summarized in Section 4.2.4. 
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Table 10. Approximate analytical expressions for the neutrino mass observables for the extreme 
cases of the mass ordering. For Ovfif) the typical (isotope-dependent) half-lifes are also given. 





E 














NH 


~ 0.05 cV 


^sin^eiaAml + IC/esPAm^^ 
~ 0.01 cV 


sin2 0i2^Am| + |(7e3iyAm2^e2'(— « 
~ 0.003 eV => T^J^ ^ 1028-29 yj.g 


IH 


~ 0.1 eV 


^Ami 
~ 0.05 eV 


^ Amj^ \/l — sin^ 29i2 sirP a 
~ 0.03 eV T0^2 ^ 1026-27 


QD 


3mo 


mo 


moy^l — sin^ 29i2 sin^ a 
> 0.1 eV => TO^^j > 1025-26 



Table 11. "Neutrino mass matrix" for the present decade. It is assumed that KATRIN will reach its sensitivity limit of = 0.2 
eV, that Oi//9/3-experiments can obtain values down to {niee) = 0.02 eV, and that cosmology can probe the sum of masses down to 
E = 0.1 eV. N-SI denotes non-standard interpretation of Oi^/9/3, N-SC is non-standard cosmology. 



KATRIN 0^/3(3 cosmology 







yes 


no 


yes 


no 


yes 


no 


KATRIN 


yes 






QD + Majorana 


QD -1- Dirac 


QD 


N-SC 


no 






N-SI 


low IH or NH or Dirac 


m„ ^ 0.1 eV or N-SC 


NH 


0ul3l3 


yes 


• 


• 






(IH or QD) -1- Majorana 


N-SC or N-SI 


no 


• 


• 






low IH or (QD + Dirac) 


NH 


cosmology 


yes 


• 


• 


• 


• 






no 


• 


• 


• 


• 







The transition to the quasi-degenerate regime takes place when 7713 ^ 0.03 eV. 
If the smallest mass assumes such values, the normal and inverted mass order- 
ing generate identical predictions for the effective mass. The results in this case 
are therefore identical to the ones for the normal mass ordering treated above in 
Section 4.2.1. 

For the inverted hierarchy case mp is again below KATRIN's sensitivity, and 
S ~ 2y^ Am\ ~ 0.1 cV is in the range of future limits. Table 10 shows the mass 
observables for the IH scheme. The values of the effective mass for the various 
special cases are displayed in Fig. 13. In Table 11 it is attempted to illustrate 
the complementarity of neutrino mass observables. Prospective sensitivity values of 
mp = 0.2 eV, (mee) = 0.02 eV, and S = 0.1 eV are assumed and the interpretation 
of positive and/or negative results in all 3 approaches is given. 

4.2.3. Mass scale 

As mentioned above, from the fundamental quantity (mee) one can also extract 
information on the masses of the individual neutrinos. We focus here on the quasi- 
degenerate regime, which is the easiest, though still non-trivial, case. The smallest 
effective mass can be written as 

(-ee)Sf„ - rno (iC/etP - |C/e.P ~ |f/e3p) = -o ^ ' ^^^^ „ ^ | t/e3 P ^ ^^^^ 
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Fig. 14. Lines of constant mo in the sin 612 — sin 613 plane, predicted for a QD mass spectrum 
with (rriee) = 0.1 eV. Also shown is an allowed 3(j region in sin^ 812 and sin^ diz- 



We show in Fig. 14 iso-contours^^^ of toq in the plane spanned by sin^ 6*12 and 



|[/e3p- With a hmit (wee)°^n on the effective mass at hand, one can translate this 
into a limit on the neutrino mass, which reads 

1 + tan^ 6*12 _ , 



Too < ("T-ee) 



cxp 

™" l-tan2 0i2-2|C/e3| 



\cxp 
/ min 



12, t7l3 



(52) 



This function /(6'i2, 9i^) varies from 2.57 to 3.29 at Icr and from 2.17 to 4.77 at 3(t. 
The limit on the effective mass is about 0.5 eV (see Table 2), and hence toq < 1.6 
eV and 2.4 eV, respectively. Therefore, the current limit on toq from Qvpp is very 
similar to the one from the Mainz experiment. 

Perhaps more interesting is the determination of the neutrino mass scale in 
future experiments if information from complementary neutrino mass observables 
is combined. For instance^^'^, consider the scenario defined by 
TO3 [eV] (m,,) [eV] mp [eV] S [eV] 



0.3 



0.11 - 0.30 



0.30 



0.91 



The prospective errors one can use are cr(TO^) = 0.025 eV^ and cr(S) = 0.05 
eV, and an "experimental error" cr{{mf,e) c^^) = \ {rriee) e^^cr{Tohs) /^ohs, where 
o-(robs)/robs is motivated by the GERDA proposaP^ to be ~ 23%. The "the- 
oretical error" from the NME uncertainty was defined as a{{mee)) = (1 + 
^) ^(m-ee) + cr((TOee)pxp)j — {rriee) ■ Depending on the measured effective mass 

(TOee)gxp O"^^ '^^^ ^O'^ obtain the values of toq which can be reconstructed. Fig. 15 
shows the results of the analysis. If C = one finds alm^) ~ 15% at 3a, while for 
C — 0.25 it holds that (t(to3) ~ 25%. If one includes a wrong cosmological input 
the reconstruction of 7713 can be wrong by up to one order of magnitude. Leaving 
S out of the analysis yields (t(to3) ~ 50%, showing that the precision is largely 
determined by cosmology. 

A detailed analysis was performed in Ref.^'''^, from where wc have taken Fig. 16. 
As was noted in that paper, the uncertainty of the oscillation parameters is of 
little importance in determining toq. To take into account the NME uncertainty 
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Zexp=0.3eV 
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QD; 




E m3=0.3eV 
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/m=0.3eV 
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0.1 

1113 |eV] 
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0.1 0.3 
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Fig. 15. 1(T, 2a and 3(t regions in the ra^- (mee) ^-^p plane for a quasi-degenerate neutrino mass 
scenario. The upper plots are for no NME uncertainty, the lower plots assume 25% uncertainty. The 
left plots show the correct (solid line) as well as two possible incorrect cosmological measurements 
(dashed lines). The right plots leave S out of the fit. The area denoted HDM is the range of 
(mee) from the claim of part of the Heidelberg-Moscow collaboration. Taken from^^''. 



the following procedure was proposed: A4 is the unknown true NME and A4o is the 
NME used to obtain (rriee) , which denotes the effective mass extracted from an 
experiment. The parameter F is connected to the ratio = Ai /Aio in the sense that 
^ ranges from l/VF to F. If the experimental error on {mee)axp sufficiently small 
0.06 cV for NME uncertainty F ^ 3), the neutrino mass spectrum will be shown 
to be QD, and toq will be constrained to lie in a rather narrow interval of values 
limited from below by mo ^0.1 eV. The uncertainty in the NME directly translates 
into an uncertainty in mg, in analogy to Eq. (52). In the case of an intermediate 
value of (mee)g^p = 0.04 eV, shown in the middle column of Fig. 16, an allowed 
range of 0.01 eV ^ toq ^ 0.1 eV could be established for precise measurements. In 
the case of an inverted ordering only an upper bound toq ^0.1 eV will be obtained. 
This result can be easily understood from the usual (mee) vs. smallest mass plots, 
which arc basically flat for the inverted ordering and ^0.1 eV. 
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Fig. 16. The reconstructed range for the lightest neutrino mass at 2(j C.L. for normal and inverted 
mass ordering as a function of the Icr experimental error on ('Tiee);,j,p (here called | <m> |°'"'). 
The results are shown for three representative values {mee)°^^ = 0.004, 0.04, 0.2 eV, and for fixed 
NME (first row), and an uncertainty of a factor of _F = 2 and F = 4 in the NME (second and 
third rows). The dashed (solid) lines correspond to the present uncertainties in the oscillation 
parameters. To the left of the dotted lines, a positive signal is obtained at 2a, whereas to the right 
only an upper bound can be set. Taken from^^^. 

Other analyses on neutrino mass extraction from different neutrino mass exper- 
iments including Oi^/3/3 can be found in Refs.^^^'^^'^^^. 

4.2.4. Mass ordering: testing the inverted hierarchy 

From Fig. 10 the interesting possibility of ruling out the inverted mass ordering 
becomes obvious. The minimal value of the effective mass, repeated here for con- 
venience, is non-zero and given by 

MZn = (1 - iC^eaP) v/a< (1 - 2 sin^ 9,^) . (53) 

If a limit on the effective mass below this value is obtained, the inverted ordering 
is ruled out if neutrinos are Majorana particles. In case the mass ordering is known 
to be inverted (e.g. by a long-baseline experiment or by observation of a galactic 
supernova) then the Majorana nature of neutrinos would be ruled out. 

As a rough requirement for experiments we calculate the difference between the 
minimal effective mass for the inverted ordering and the maximal effective mass for 
the normal ordering multiplied with the nuclear matrix element uncertainty factor 
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,02 sitfe|,=0.33 



0.(11 (1,02 0.03 0.04 0,05 



Fig. 17. The difference A(mee) of {'"ee)I"j',i and C('"ee)5^°ax ^ function of sin^ S13 for differ- 
ent nuclear matrix element uncertainty factors f = 1,1.5,2 and 3 (from top to bottom). We 
have chosen an illustrative value of the smallest mass of 0.005 eV and sin^ 612 = 0.28 (left) and 
sin2ei2 = 0.33 (right). 



►160,152,153 



^ mv 
' mill 



C (mee) 



nor 
max ■ 



(54) 



We plot this difference as a function of sin^ ^^13 in Fig. 17. Obviously, the largest 
dependence stems from 9i2, and the smaller 612 is, the better. This is clear from the 
previous discussion and Fig. 13, because the smaller 9i2 is, the larger is {nT-eeY^^- 
The effect^^"^ of non-zero 613 is to shghtly decrease (mee)™i'n ^^'^ *o shghtly increase 

/ \ nor 
\TOee/iii ax- 
One can translate the effective mass necessary to rule out (or touch) the inverted 

hierarchy into half-lifes. The very important dependence on 612 has recently been 

discussed in Ref.^''. The plots in Fig. 18 are generated using the compilation of 

NMEs from Table 5. The current 3(t range corresponds to an uncertainty of a 

factor 2 in the minimal value of the effective mass, which is of the same order 

as the current uncertainty in the NMEs. The factor 2 due to 612 corresponds to 

a factor of 2^ = 4 in half-life. In experiments with background, see Eq. (6), this 

means a rather non-trivial combined factor of 2'^ = 16 in the product of measuring 

time, energy resolution, background index and detector mass. Therefore, a precision 

determination of the solar neutrino mixing angle would be very desirable to evaluate 

the requirements and physics potential of upcoming Oi//3/3-experimcnts in order to 

test the inverted ordering^^. 



4.2.5. Majorana CP phases 

Apart from measuring the effective mass in case of a normal hierarchy, determina- 
tion of a Majorana CP phase from ncutrino-lcss double beta decay is probably the 
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Fig. 18. Required half-life sensitivities to exclude and touch the inverted hierarchy for different 
values of 612. The upper plots show the necessary half-lifes for sin^ 9i2 = 0.27 (upper left) and 
sin'^ 612 = 0.38 (upper right). The lower left plot includes the current 3cr uncertainty for 612. 
The lower right plot shows the necessary half-lifes in order to touch the inverted ordering, which 
is independent on 612. The small horizontal lines show expected half-life sensitivities at 90% 
C.L. of running and planned Oi^/3/3-experiments. When two sensitivity expectations are given for 
one experiment they correspond to near and far time goals, see Table 4. Taken from^S. 



most difficult physics goal related to this process'. One general point to be made 
here is that there is only one observable, {rriee), and thus only one of the two Ma- 
jorana phases (or a combination of the two phases) can be extracted. In addition, 
complementary information on the neutrino mass scale has to be put in for such a 
measurement. A final remark is that the process is not CP violating, i.e. the rate 
of the Qv^^ [3^ process depends on the same quantity as the Oi//3/3 process-*. 

The prospects of measuring the CP phase in neutrino-less double beta de- 
cay have been discussed in several papers^^'^. A somewhat pessimistic conclusion 



'Though |f/e3| has some influence on f)ul3j3^^^, extracting it from a measurement is also 
inpractical^^^ . 

■■Manifest CP violation from Majorana phases is discussed e.g. in^^^. 
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sin^Bj^ = 0.25 ± 3% sin^Gj^ = 0.31 ± 3% sin^Oi , = 0.38 ± 3% 




uncertainty in l<m>l from NME 



Y///A data consistent with = 7i data consistent with a.,j = 

I \<m>\ and Y. inconsistent at 2a CP violation established at 2a 

sin^Ojj = ± 0.002, Am^j = 8x10'^ ± 2%, Anjjj = 2.2x10"^ ± 3% observed l<m>l = 0.3 eV 

Fig. 19. Constraints on the Majorana phase la (here called 021) at 95% C.L. from an observed 
{mee)(,xp = 0-3 eV and prospective data on S, as a function of the NME uncertainty factor F. 
Shown are the regions in which the data are consistent with one of the CP-conserving values 
(hatched), observed S is inconsistent with ('Tiee)j,^p (light-shaded), and Majorana CP-violation is 
established (red/dark-shaded). Taken from^^^. 

has been drawn in Rcf.-'^^^, whereas the requirements for such a measurement have 
been discussed in^^^, and found to be not too unrealistic. 

The requirement for determining the phases is clear from Figs. 10 and 12. Ex- 
perimentally one should find results lying in the areas indicated with "CPV" , which 
are however smeared by experimental and theoretical uncertainties. This is realistic 
only for the inverted ordering or the quasi-degenerate scheme. Neglecting ^13, the 
effective mass is in these cases is proportional to 

(mee) oc jcos^ 6*12 + 6^*" sin^ 6*121 = \j 1 - sin^ Id^i sin^ a . (55) 

Therefore, the larger Q\2 is, the more promising it is to extract a from measure- 
ments. Recall that ruling out the inverted mass ordering is easier if ^12 is small. 

A detailed statistical analysis has been performed in^^-'^, from which we present 
Fig. 19. One can see that, as expected, for larger values of 6*12 the areas in pa- 
rameter space become larger. For instance, if sin^ Qy^ ^ 0.3 and ~ 10% errors in 
the measured (''T^ee)oxp ^'^'^ ^ present, the NME has to been known to better 
than within a factor of 1.5. For smaller values of the errors, ~ 0.01 eV and 
(Ts — 0.05 eV, Majorana CP-violation could be established even for F ~ 2 See 
Section 4.2.3 for the definition of F). Finally, the Majorana phase 2a has to have 
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a value approximately in the interval ^ (7r/4 — 37r/4). In the inverted hierarchy 
the required errors have to be smaller, and the determination of the phase is more 
challenging. 

4.2.6. Vanishing effective mass 

Unfortunately, the normal mass hierarchy can allow for complete cancellation of 
the effective mass (see e.g.^^^). In terms of Fig. 8, this "cancellation regime" means 
that a triangle can be formed. If \Ue3 \ = then the requirement is 

nil „ 1 , ^ 

— = tan^ 012 ~ - , 56 
m2 2 

while for mi = one needs 

m2 tan^ Oi^ o ^ 2q /r^N 
— = . o ^ ^ 3 tan 6li3 . (57) 
ma siir 6*12 

In both cases the Majorana phases need to be such that the two surviving terms 
have opposite sign. For the case of arbitrary 6*13 one finds^^^ 



cos 2a 



^3 '^13 C]^3(mj^ C]^2 + '^2 ^12) 



2 mi m2 s?2 c?2 cf3 ' /^gN 

sf^ + 4^ {ml sj^ - ml 

COo ZijD — 9 9 9 

2 m2 ma sjj sfg cf^ 

It may seem unnatural that the 7 parameters on which {m^^) depends conspire in 
such a way that the effective mass vanishes. However, recall that the effective mass 
is the ee element of the Majorana neutrino mass matrix'*. This matrix is generated 
by the underlying theory of mass generation, and texture zeros occur frequently 
in such (flavor) models, see Ref.^''° for a general analysis and^^^ for symmetries 
leading to (mee) = 0. 

One may ask whether the effective mass remains zero, or whether corrections 
lead to small but non-zero (mee). In fact, there are several possibilities for non-zero 
0://3/3-rates, even if (mee) = 0: 

• the first point to make here is that the dependence of the amplitude on the 
neutrino parameters goes as (see Eq. (40)) 

While the second term is very much suppressed by mf/q^ ^ 10"^^ with respect 
to the usual effective mass term, it is in general non-zero, even when (mee) is 
zero; 

• another source of correction arises from radiative corrections. While in the 
effective theory the rcnormalization of the mass matrix is multiplicative (see 

''Zeros of the remaining elements of have been studied in^^^, the presence of two zeros in 
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Section 4.2.7), this may no longer be true in the theory beyond the effective 
one. For instance, if rtii, is generated via the type I seesaw mechanism, then 
threshold corrections from integrating out the heavy neutrinos one by one will 
in general lead to non-zero {niee), even if the ee entry of mj^M^^ niD is zero; 

• in the type I seesaw mechanism, mjj/Mn is actually only the leading order 
term. "Ncxt-to-leading order" corrections m'^/M'^ are present and in general 
induce non-zero terms in the ee entry, even if {mj^ -^r^ fnD)ee is zero^^^; 

• if another source of lepton number violation is present, then a non-zero ee entry 
of the neutrino mass matrix will be induced via the Schechter-Valle diagrams 
from Fig. 1, even if {rriee) ~ 0; 

• finally, it is plausible that a flavor-blind Planck scale term is present, which 
induces an effective mass of order u^/Afpi ~ 10~^ eV. This term arises from 
the Weinberg operator Eq. (17) with the Planck scale inserted as A. 

All these sources give of course very small but in general non-zero contributions 
to the effective mass. One might ask whether one can determine experimentally by 
other means if the effective mass vanishes. While this is not possible, one can show 
however that the effective mass cannot vanish: from Fig. 12 note that (TOee) — 
corresponds to to^ ^ 0.02 cV and S ^ 0.1 eV. Thus, finding these quantities above 
such values immediately rules out the possibility of vanishing (mee)- Of course, 
determining experimentally that the inverted mass ordering is realized also implies 
that (mee) 7^ 0. 



4.2.7. Renormalization 

The renormalization group (RG) evolution of neutrino parameters has recently been 
reviewed in^'''^. If some unknown high energy theory at a scale A leads to a mass 
matrix m^, then in the effective theory one has the following mass matrix at low 
scale A, where measurements take place: 

■(m°)ee/e (TO°)ep/e/^ 

Ic. I • {ml)^rl^.lr 1 , (59) 




where 



C A 1 A 

lOTT^ A lOTT^ A 



with a,P G {e, fi, r}, C = 1 in the MSSM and C = -| in the SM. One can safely 
drop Ue and from the above expression and describe the RG evolution with I^- 
and /q,^ only. We furthermore have 

= -3ff| + 2(y? + yf, + y^) + 6 (y? + y^ + yl + yl + yl + yl) + Ah , . 

^MSSM ^ _y2 _ + 6 (j,2 + y2 + ^2) _ (blj 

Here 31^2 are the electroweak gauge couplings, yx the Yukawa coupling of fermion x, 
and Ah the Higgs self-coupling. The RG evolution of {mee) is therefore basically a 
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Fig. 20. Extrapolation of the effective mass from 0.35 eV at low scale to fiighcr energies. Tfie SM 
curves in the left plot correspond to Higgs masses of 114 GeV, 165 GeV and 190 GeV (from bottom 
to top). In the MSSM plot on the right, the Higgs mass is 120 GeV, tan ^ = 50, MsuSY = 1-5 TeV. 
Taken from^'^"'. 



rescaling of the effective mass with la,^ . In contrast to the running of the individual 
parameters of (mee) (^12, ^13, ttii, 7712, TO3, a and /3), which can be very dramatic, 
the RG evolution of {m^f.) is modest. Its running does basicahy not depend on the 
mass ordering or any of the other neutrino mass and mixing parameters. It is an 
interesting exercise to consider the j3 functions of the 7 parameters of {m^f,) and 
to show that at the end all dependence on ^12, ^13, mi, 7712, ms, a and /3 drops 
out. The effective mass typically increases from low to high scale. Fig. 20 shows an 
example for its running^'''*. 



4.2.8. Distinguishing neutrino models 

We have mentioned in Section 4.1.2 that the peculiar and unexpected form of 
lepton mixing (see Eq. (31)) is assumed to have its origin in the presence of flavor 
symmetries^ There is a large abundance of such models, many leading to the same 
neutrino mixing scheme, for instance tri-bimaximal mixing (TBM). The question 
arises how to distinguish them from one another. It turns out that neutrino mass 
observables can help in disentangling the vast amount of flavor symmetry models. 
One example is that the flavor symmetry leads to correlations of the mass matrix 
elements, which imply correlations of observables. For instance, the effective mass 
could be correlated with the atmospheric neutrino parameter sin^023, which was 
obtained in a model in Ref.^''^, see Fig. 21. Recall that in general ^23 has no influence 

on (mee). 

Another point are "sum-rules" : the most general neutrino mass matrix giving 
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Fig. 21. Left: correlation between the effective mass and the atmospheric neutrino mixing angle in 
a specific flavor symmetry model. Taken from^^^. Right: allowed regions in {rriee) — E parameter 
space for the sum-rules — H — — = — (top) and — — I — — = — (bottom), for both the TBM 

^ m2 m3 mi ^ ^' mi 7712 m^ ^ 

(black) and 3cr values (light red) of the oscillation data, as well as for the sum-rules violated by 
30% (green hatches). Taken from^'^^. 



rise to TBM is 



B 

-B 



V 



B 

-B 
-B 



D) 



(62) 



As such, the (eomplex) eigenvalues A — B, A + 2 B and D are independent of the 
mixing angles: no matter what A, B, D are, the PMNS mixing is given as Eq. (31). 
However, very often the structure of the mass matrix is simpler than in Eq. (62), and 
"sum-rules" between the neutrino masses arise. Examples are^^*" 2/to2 + l/ms = 
1/mi or 1/to2 + 1/^3 = 1/wi, and detailed studies of the predictions can be 
found in Ref.^^^, from which we took the right plot in Fig. 21. Other discussions 
on mass-related phenomenology of flavor symmetry models can be found in^^^. 



4.2.9. Light sterile neutrinos 

The easiest way to depart from the standard 3 neutrino picture discussed so far is 
to add light sterile neutrinos. In fact, we have mentioned in Section 4.1 several hints 
which point to the existence of additional radiation in the Universe, as well as for 
one or two additional mass-squared differences in the eV regime. With one or two 
light sterile neutrinos the PMNS matrix becomes a unitary 4 x 4 or 5 x 5 matrix. 
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Short-baseline oscillations depend on UeA (and Ueb) as well as (and Only 
eV-like mass-squared differences play a role for such experiments, and only if two 
sterile states are added does the possibility of CP violation arise. This can explain 
the different neutrino and anti-neutrino results from MiniBooNE and MiniBooNE 
plus LSND, respectively. 

Table 12 shows the results-'^^® from a global fit to the world's short-baseline data, 
taking into account the recent re-evaluation of reactor fluxes^^^. The data are not 
sensitive to whether the two sterile neutrinos are above or below the three active 
ones (2+3 or 3-1-2 scenarios), but are sensitive to whether the active neutrinos are 
sandwiched between two sterile ones (1+34-1). Note that the fit to 1+3+1 scenarios 
is slightly better than for 3+2/2+3 scenarios. Fitting only the reactor experiments 
is possible in a 3+1 or 1+3 scenario, and gives^^* \Uei\ = 0.151 and /S-mli = 1.78 
eV2. 

If there are two sterile neutrinos, the nomenclature for the 8 possible mass 
orderings is as follows: 

(i) SSX, where X = N for a normal and X = I for an inverted ordering of the mostly 
active neutrinos. In these schemes the two predominantly sterile neutrinos are 
heavier than the three predominantly active neutrinos (2+3 scenarios); 

(ii) XSS (X = N or I as before), where the two predominantly sterile neutrinos are 
lighter than the three predominantly active neutrinos (3+2 scenarios); 

(iii) SXS with X = N or I, where the three active neutrinos are sandwiched between 
the sterile ones (1+3+1 scenarios). In this class there can be four possible 
scenarios, which we denote as SXSa and SXSb. The scheme SXSa corresponds 
to the state 1^5 higher than the three active states and SXSb corresponds to the 
state vc, lower than the three active states. 

The individual masses, expressed in terms of smallest mass and the four mass- 
squared differences, can be found in Ref.^^^. The new sterile neutrinos will con- 
tribute to the sum of masses in cosmology, to the kinematic mass in KATRIN, and, 
if they are Majorana particles, to Oz^/3/3. The effects of sterile neutrinos on neutrino- 
less double beta decay have been studied by various authors^^^'^^^'^*^'^^^'^^". One 
simply extends the sums in the definitions of E, and (mee) from i = 3 to i = 5. 
The interpretation of reactor experiments actually observing oscillations, which is 
possible after the new reactor fluxcs^^^ are taken into account, makes the applica- 
tion of the results to Qvfifi easier. If only LSND and MiniBooNE supply oscillation 

data, the 1^^^— >-W transitions depend on Uei and J7^i {i = 4,5). For mp and (mee) 
only Uei is required, and one needs to assume something about U^i to extract Uei 
from the fit results. However, reactor oscillation survival probabilities depend only 
on Uei- There are two more Majorana phases which show up in the modified effec- 
tive mass, which is the sum of the contribution considered so far ((mge)^'^') plus 
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Table 12. Parameter values and x at the global best-fit points for 3+2 and 1+3+1 
oscillations. Taken from^^®. 

Am|,[eV2] |j7e4| |C/^4| Am|,[eV2] {U.^l l^/^sl xVdof 

3+2/2+3 0.47 0.128 0.165 0.87 0.138 0.148 110.1/130 

1+3+1 0.47 0.129 0.154 0.87 0.142 0.163 106.1/130 



new terms from the sterile states ((mee 



The usual phase space factors and matrix elements as in the standard interpretation 
apply (nuclear physics is the same as for the standard case, as long as the masses 
do not exceed q ~ 100 MeV). The additional contribution {meeT^ from the sterile 
neutrinos could be leading, sub-leading, or of the same order of magnitude as the 
active neutrino part (mee)^"^*- 

Neglecting the smallest neutrino mass and using the best-fit values from Table 
12 gives the following predictions for the mass observables^*°: 

SSN: the active neutrinos give the same contribution as in the standard case for 
NH. The contribution to the sum of masses from the sterile states dominates 
and is given by E = Am|j^ + ^JKm^ ~ 1.62 eV. The contribution to mp 
is ■\/ATO4]^|[/e4p + Amgj^|[/e5p — 0.16 eV. The contribution to the effective 
massis \UeA? ^AmJ^ + e^'^*!"*^' \Ueb? , which is between 0.007 and 

0.029 eV, and hence larger than the typical value for NH. Thus, the effective 
mass cannot vanish (for the best-fit point), in contrast to the standard case. 
SSI: the active neutrinos give the same contribution as in the standard case for IH, 
and the sterile states give the same predictions for the mass observables as in 
SSN. The effective mass can therefore vanish, in contrast to the standard three 
neutrino case. 

NSS: the active neutrinos are QD (normal ordering) with a mass scale ^jKm^ ~ 
0.93 eV, and this governs the predictions for {rriee) and mp. For cosmology, 
S = 3 y/Aml^ + y/Aml^ - Aml^ ~ 3.4 eV. 
ISS: same as NSS, except for inversely ordered active neutrinos. 



SNSa: the active neutrinos are QD (normal ordering) with a mass scale a/ Am\^ ~ 0.69 
eV, and this defines the predictions for {rriee) and mp. The sum of masses is 

E = 3 a/Ato|i + y/Aml^ + Amj^ ~ 3.2 eV. 

SNSb: the active neutrinos are QD (normal ordering) with a mass scale ^/Amf^ — 0.93 
eV, and this defines the predictions for {niee) and mp, up to a small correction 
of order |J7e4p \/ Arri^^ + Am|j^ ~ 0.03 eV e.g. for {rriee)- The sum of masses is 
E = 3 y/Amj^ + y/Amj^ + Amj^ ~ 4.0 eV. 
SISa: same as SNSa, except for inversely ordered active neutrinos. 
SISb: same as SNSb, except for inversely ordered active neutrinos. 
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Note that these are all lower limits, because we have put the smallest neutrino 
mass to zero. In case of SSI the effective mass can vanish when the 3 active neutrinos 
are inversely ordered, in contrast to the three neutrino case in Section 4.2.2. If 
future Oz^/3/3-experiments measure a tiny effective mass, or obtain a limit below 
(mee)min giveu in Eq. (50), and the neutrino mass ordering is confirmed to be 
inverted from long-baseline neutrino oscillations, the sterile neutrino hypothesis 
would be an attractive explanation for this inconsistency. This is the first example 
in which a deviation from the standard picture of 3 active neutrinos shows up and 
influences the interpretation of Oz^/3/3. We show in Fig. 22 the effective mass against 
the smallest mass for the 1+3 and 2+3 cases. 

Obviously all schemes have difficulties with cosmology, the contribution to the 
sum of masses exceeds 1.5 eV in all cases. KATRIN and next generation QvjijS- 
experiments will see a signal in all cases except for SSI and SSN, unless the masses 
and mixings take values at the very high end of their currently allowed ranges. An 
analysis of KATRIN's potential to separate one or more sterile neutrino component 
from the active neutrino component has been performed in-'^^^. It was shown that 
KATRIN will definitely be able to separate one or more sterile neutrino components 
from the active neutrino ones, if they do in fact have mass and mixing in the range 
considered here. With a limit on the effective mass being around 0.5 eV (see Table 
2), the schemes with QD neutrinos with mass scale yAmf^ or \/Km^ are in fact 
already tested, giving constraints on the Majorana phases a already at the current 
stage^^^. 

If the "reactor only" results of |?7e4| = 0.151 and Am|j^ = 1.78 eV^ are used, 
then E > Ato|j^ ~ 1.3 eV or E > 3 \/ /S.m\^ ~ 4.0 eV, depending on whether a 1 
+ 3 or 3 + 1 scheme is realized. The contribution to KATRIN is either 0.52 or 1.3 
eV, and the effective mass either receives a contribution of 0.03 eV, or corresponds 
to QD neutrinos with a mass scale y'Krn^. For the 1+3 case, again, the effective 
mass can vanish if the active neutrinos are inversely ordered. 

The existence of sterile neutrinos can also be tested in upcoming oscillation 
experiments and via cosmological observations, as discussed in Section 4.1.2. 

4.2.10. Exotic modifications of the three neutrino picture 

Exotic modifications of the 3-neutrino framework are possible and may spoil the 
discussion presented so far in this Section. 

The most obvious modification is that neutrinos are Dirac particles, in which 
case there is no neutrino-less double beta decay and writing this review was all 
in vain. A useful way to show this in the effective mass is to note that one Dirac 
neutrino can be written as two maximally mixed Majorana neutrinos with common 
mass rrii and opposite CP parity. The effective mass is then 




(63) 
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Fig. 22. Top: Effective mass against the smallest mass for the 1+3 scheme, with Am|j^ = 1.78 
eV2, |[/e4p = 0.023 (dark) and the 2a range Am|^ = (1.61 - 2.01) eV^, \U^i\^ = 0.006 - 0.040 
(light). Bottom: same as above for 2+3 scheme, with Am^-^ = OAT eV^, Arrigj^ = 0.87 eV^, 
\Ue4\^ = 0.016, |(7e5p = 0.019 (dark) and Am|^ = (0.42 - 0.52) eV^, Am^^ = (0.77 - 0.97) eV^, 
|C/e4|^ = 0.004-0.029, |(7e5p = 0.005-0.033 (light). The black solid and dashed lines correspond 
to the standard 3 neutrino best-fit and 2a cases. 



A small splitting of the degeneracy can be described with the mass matrix 

+ and ^?=^. (±1 + 1), (64) 

with the indicated new eigenstates and mixing matrix. These Pseudo-Dirac neu- 
trinos lead to a contribution to the effective mass of about enii = ^ 5m^ /rrii, with 
Sm^ = [mfY ~ {^iY- Regarding limits on such splitting, roughly speaking, values 
larger than 5m? ~ 10~^^ eV^ for mi and m2 are forbidden by solar neutrino data, 
and Sm^ ^ 10^"^ cV^ for 7713 by atmospheric data^^^. If all three states are Pseudo- 
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Dirac the effective mass is basically zero"'^*^, while if one or two are Pseudo-Dirac 
interesting predictions for the effective mass arise. This can happen in "bimodal" 
or "schizophrenic" scenarios^^^, in which at leading order one or two mass states 
are Dirac particles while the other one is Majorana. Because lepton number is not 
conserved, loop corrections imply small Pseudo-Dirac terms for the Dirac states. 
For instance, if V2 is a Dirac particle then the effective mass in the inverted hierar- 
chy isis-* (mee) ^ \/Ami cf 2 cfg, roughly a factor of two larger than the minimal 
value in the standard case, see Eq. (49). A generalization to all possibilities can be 
found in^®'^. 

A possible modification of the three neutrino picture mentioned before is the 
possible non-unitarity of the PMNS matrix, which has however negligible effect 
on the effective mass^^^. 

Another exotic property is CPT violation. Interesting consequences for 
Oz//3/3 have been considered in^^^, where a simple one family example is discussed. 
In the (z^, v) basis, where CPT transforms v into i/ up to a phase, the mass matrix 
can be written as 

M=^''t^/^v (65) 



y M - A ^ 

Here y mixes v and v, while A leads to different masses for v and v. The eigenstates 
v± with masses m± = /i ± -^/jj/p + can be shown to be Majorana neutrinos 
(i.e. CPT transforms v± into v± up to a phase) only if A = 0. This in turn would 
imply however that CPT is conserved. CPT is violated for A 7^ 0, in which case 
neutrinos cannot be Majorana particles. The amplitude for Qvf5fi sums over 771+ 
and m_ and is non-zero Therefore, neutrino- less double beta decay takes place 
even if neutrinos are strictly speaking not Majorana particles. The neutrino-less 
double positron decay proceeds with the same "effective mass" . 

In principle can also provide limits on parameters associated with violation 
of Lorentz invariance or the equivalence principle. The constraints^^^ on the 
difference of maximal velocities of mass states or on non-universal couplings of 
neutrinos to the gravitational potential arc in general weaker than the ones from 
neutrino oscillations^^'^. 

It should be mentioned here that 2i^/3/3 constrains violation of the spin- 
statistics theorem for neutrinos. With two identical particles in the final state 
there are two diagrams with exchanged momenta pi -H- j32 • Their relative sign de- 
pends on whether Fernii-Dirac or Bose-Einstein statistics applies. By writing the 
amplitude A as cos^ -^fomiionic + sin^ ^-^bosonic, conservative limits of sin^ (p ^ 0.5 
can be set^^^. 



5. Non-Standard Interpretations 

After discussing is some detail the standard interpretation of neutrino-Icss double 
beta decay, we turn to non-standard interpretations, repeated here for convenience: 
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Neutrino-less double beta decay is mediated by some other lepton number vi- 
olating physics, and light massive Majorana neutrinos (the ones which oscillate) 
potentially leading to OvPfi give negligible or no contribution. 

It is convenient to express the decay width of neutrino-less double beta decay in 
the following form (see Eq. (4)) 

r""^ =^G,(Q,Z) |7W,77,^ (66) 

X 

Here we sum over all possible mechanisms which arc denoted by a subscript x, 
with matrix element Mx and a dimcnsionless particle physics parameter rjx- For 
the standard interpretation of light neutrino exchange, 

?7i = (mee)/me <9.9 X 10-^ (67) 

Note that different mechanisms can interfere coherently, a case we will discuss 
in Section 6.4. Most of the times the alternative mechanism is connected to a 
high energy scale. The corresponding particle physics amplitude, which has to be 
compared with the standard one G\ {m^^ jq^ from Eq. (40), could be written as 

v4hoavy ^ • (68) 

Here c contains new Yukawa and/or gauge couplings and A is the new physics 
scale. This is a helpful but crude approximation, which is in fact not fulfilled by 
several mechanisms to be discussed in the following. However, the current limit 
(rriee) — 0.5 eV corresponds to A ~ TeV, by all means an interesting energy scale. 
In fact, we will encounter in what follows some alternative mechanisms with LHC 
phenomenology. On the other hand, it means that if the new physics scale exceeds, 
say, 10 TcV, then it will not contribute significantly to 0i//3(i. In what follows we 
will aim at a complete list of non-standard realizations of neutrino-less double beta 
decay, for earlier reviews see^^^'"'^^^. 

An ideal experimental signature for drawing the conclusion that a mechanism 
different from active neutrino exchange is present would be that KATRIN and 
cosmological observations do not see a signal, but Oi^/3/3 is observed with a half-life 
corresponding to, say, (mge) — 0.5 eV. To put in another way, in plots of neutrino 
mass observables, such as in Fig. 12, one ends at points outside the allowed areas. 
On the other hand, if one ends in those plots in the allowed areas, then it is not 
necessary to consider non-standard interpretations, except for setting limits on the 
associated parameters. 

5.1. Heavy neutrinos 

An interesting way of realizing Oj^/3/3 is through the exchange of heavy Majorana 
neutrinos^^^'^^'^'^^^. The Feynman diagram is the same as in Fig. 8, with the neutri- 
nos not being the ones whose oscillations are observed, and with the PMNS matrix 
elements Uei replaced by Sei, where S is the matrix describing the mixing of the 
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rate 



i 






" mass 



Fig. 23. Left: "lobster" diagram for lepton number violating processes with Majorana neutrino 
exchange. Right: typical behavior of Of/S/S-like processes with free Majorana mass m: for small 
masses m? <^ the rate increases with m^, while for large masses m? ^ it decreases with 
. Here "small" and "large" are defined relative to the energy scale q^ of the process, which 
could be the mass of a decaying particle/nucleus, or the center of mass energy of a collider process. 
The maximum rate can be expected when the mass corresponds to the available energy. 

heavy neutrinos with the SM charged leptons in the charged current term. Recah 
the form of the Oz^/3/3-amplitude on the particle physics leveh 



As mentioned before, due to the typical structure of the diagram, symbolically 
displayed in Fig. 23, one sometimes calls it "lobster diagram". With Oi^/3/3 being 
a t (and u-) channel process there is no resonance. In Fig. 23 we show the typical 
behavior of Oz^/3/3-likc processes as a function of the Majorana mass. Let us stress 
that the maximum rate can be expected when the mass corresponds to the available 
energy, i.e. about 100 MeV for Qvf5(i. In analogous processes of neutrino-less double 
beta decay (see Section 7) the energy scale and therefore the range in which the 
strongest limits on the mass and mixing arise, may be different. In addition, there 
could be s-channel processes in which a resonance could be hit, leading to even 
stronger constraints. If <; 100 MeV, the rate is proportional to 



We will focus on this case of heavy neutrinos. 

Turning to nuclear physics, the neutrino potential in Eq. (13) is modified be- 
cause the neutrino energy and momentum are dominated by its heavy mass. A 
dependence on the axial mass Ma — 0.9 GeV, which appears in the nuclear form 
factors and which takes into account the finite size of the nucleons. is introduced 




(69) 




(70) 
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because of the short-range nature of the process. Without the form factor the pro- 
cess would be exponentiahy suppressed due to the repulsion of the two decaying 
nuclei. Note that this is the first diagram for which is purely short-range or 

point-like. Details of the nuclear physics can be found in Refs.^^^'^'^^'^^®'^^^. Often 
one writes the contribution of heavy neutrinos as vT''^ oc {—) M\ F{A,mi)^ where 



F{A,mi) = 0(0.1) is a mildly varying function. One can write the decay width as 
= GhiQ, Z) |77h A^hP, where in the context of Eq. (66) one can write the LNV 
parameter for heavy neutrino exchange as nip {^), with nip the proton mass. The 
same phase space factors as in the standard ease apply, and the matrix elements 
absorb now various factors such as the dependence on AIa or F{A,nii). The parti- 
cle physics parameter (^) contains all singlet fermions coupling with SM charged 
leptons in the charged current terms, including heavy neutrinos from the type I and 
III seesaw mechanisms, as well as generalizations thereof, such as inverse seesaw^^^. 

Ref.^^ has recently calculated within the QRPA approach the NMEs J\4h in 
the above convention and found a range of roughly a factor of two: 172 - 412 for 
^^Ge, 165 - 408 for ^^Sc, 185 - 404 for ^°°Mo and 171 - 384 for "OTc. The spread 
originates from variation of gA, the nueleon-nucleon potential and the model space 
size. The NMEs seem to be much larger than the ones for the standard case, but 
as mentioned above they absorb several parameters. With the current limits on the 
half-lifes from Table 2 and the phase space factors from Table 1, we find 




'8 GeV"^ for^^Ge, 
« GeV"^ for^^Se, 
X 10-« GeV"^ for^""Mo, 



(71) 



* GeV"^ for^^"Te, 
The dimensionless LNV parameter for heavy neutrino exchange is 

77h - nip (-) < 1.7 X 10-« , (72) 
m 

where nip is the proton mass. Interestingly the best limit (^) < 1.8 x 10^^ GeV^^ 
stems jointly from ''^Ge and -'^^^Te. For heavy neutrinos the limit J2i \Sei\'^ < 0.0052 
from global fits applies, and this constraint on jS'eip is stronger for > 2.9 x 10^ 
GeV. Naively, one can simply compare the particle physics amplitudes G|n {niee) /(f' 
and G| {^). With {m^e) ^ 0.5 eV and q ~ 100 MeV it follows that {:^) <5x 10"* 
GeV~^, which is basically the same number as Eq. (71), given the NME uncertainty. 
As we will see in the following, the comparison of an alternative mechanism of 
OvPP to the standard mechanism on the amplitude level is remarkably successful, 
and gives constraints which are consistent with literature values taking the onerous 
nuclear physics aspects into account. Formulated provocatively, matrix elements are 
order one numbers with a corresponding uncertainty, and comparing the particle 
physics amplitudes introduces an order one factor uncertainty, which often is good 
enough to understand the particle physics implications of 01/(3(3. 

Fig. 24 shows the exclusion limits on mass and mixing of heavy sterile neu- 
trinos from Ref.^^^. The calculation covers all masses from keV to 10^^ GeV. As 
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expected, the limits are strongest when the neutrino mass corresponds to the en- 
ergy scale q ~ 100 MeV. 

There is an obvious source of heavy neutrinos, namely the ones from the type 
I seesaw mechanism, Nui^ with masses Mi. Note that these particles provide two 
sources of QvPj3: a direct one realized by exchange of Nfn and an indirect one by light 
neutrino exchange. However, the Nm are typically very heavy and have suppressed 
mixing S ~ mo /Mr ~ m^/mD — \J m^/Mn, therefore they lead to basically 
vanishing (^). Without any strong, instable and fine-tuned cancellations^^^ , the 
direct contribution from {ruee) is larger in seesaw scenarios^™'^^^'^^^. Within type 
I seesaw there is an exact relation 

Y^NliiM + SlM, = Q, (73) 

i 

where | X]-^ei™il the effective mass {rriee) in type I seesaw scenarios in which 
the PMNS matrix is strictly speaking not unitary and thus denoted here by N . The 
zero on the rhs of the above equation is nothing but the upper left zero in the full 
seesaw mass matrix in Eq. (18). Therefore, with Eq. (73), the limit on (mee) ^ 0.5 
eV directly translates to | ^ Sl^ Mi\ ^ 0.5 eV, which in the absence of cancellations 
is much more stringent than | ^ Sl^/Mi] ^ 1.8 x 10^^ GcV^^. If a low scale seesaw 
mechanism is applied and both the and the Mi are below 100 MeV, then there 
will be no neutrino-less double beta decay because^°^ of the exact seesaw relation 
Eq. (73). 

In type III seesaw scenarios the neutral component of the triplet plays the role 




Fig. 24. Exclusion plot in the \Uf,h\'^-'mfi plane, where and \U^h\ ^^'^ heavy neutrino masses 
and their mixing with the SM electron doublet. The shaded regions are excluded by Oi'/^/J-decay, 
by Big Bang nucleosynthesis and by SN1987A neutrino observations. Taken from^^^. 
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Fig. 25. Quark level Feynman diagram for the triplet realization of neutrino-less double beta decay. 



of the heavy neutrino in the type I seesaw and the discussion is analogous. 
5.2. Higgs triplets 

Higgs triplets contain a doubly charged scalar and can directly couple to two elec- 
trons and to two M^-bosons, giving rise to the quark level Feynman diagram shown 
in Fig. 25, first noted in'^*^^. This is the first diagram for Oi^^/3 which does not 
contain a neutrino line. We will show in this Section that in the simple version 
based solely on SU{2)l x U{1)y the triplet does not play a significant role in Qv^jS. 
In left-right symmetric theories this changes, and we will deal with this class of 
theories in the next Section. 

The SU (2) L triplet can be written as 

and the neutral component receives a vev (A*') = wl/2, which induces from the 
Lagrangian £a = ^Q;3i'^a*'''2 AL^, where La are Lepton doublets of flavor a, the 
neutrino mass matrix nii, = hv^. The vev is constrained from the electroweak 
p parameter to be less than about 8 GeV, and current limits on the triplet masses 
are around 100 GeV^^^. The particle physics amplitude for Ov/SP can be read off 
from Fig. (25) as 

A ^ ri2 ^ee VL ^ ^2 ('^i')ee _ ^2 (™ee) 

TO^ TO^ 

If the triplet was responsible for neutrino mass (type II seesaw) then ft.ee = 
(m,y)ee, which is the largest possible value of ftee ^l, unless unnatural cancellations 
of different seesaw terms take place. Comparing with the standard amplitude in 
Eq. (40) we see that the rate for triplet exchange is suppressed with respect to 
the standard mechanism by at least a factor {q/ra^)'^ ^ 10~^^ and hence not of 
relevance^°^'^°^. Nuclear physics details add some additional suppression^"^. 
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Fig. 26. Drell-Yan production of Higgs triplets at the LHC with subsequent decay into like-sign 
lepton pairs. In left-right symmetric theories this process is possible as well. 

There are additional diagrams in which one or two of the ly-bosons in Fig. 25 
are replaced by singly charged scalars A~ (see e.g.^°'') but those amplitudes are 
suppressed by a factor vl/v for each A~-quark vertex and by {my//'m/^~Y for 
each A~ propagator. In principle one can evade these constraints by adding exotic 
scalars with specific hypercharge and isospin quantum numbers^*^^'^^®. 

Finally, we should mention the possibility of Higgs triplet production at the 
LHC, which is possible up to masses of about 800 GeV^°^. The relevant diagram 
is shown in Fig. 26. If their branching ratio into leptons is larger than into a W 
boson pair, then their decay can give information on the neutrino mass matrix if in 
addition the pure type II seesaw is realized. In fact, BR(A a~ f3~) (x (my)^^, 

and an alternative method to probe Majorana neutrino properties was possible, as 
studied e.g. in209'2io,207^ Note in particular that the branching ratio for decays into 
two electrons is proportional to the effective mass. The other entries of the mass 
matrix could be directly studied as well, which is not possible with other processes, 
see Section 7. In case both the triplet and the type I seesaw are at work, the exact 
seesaw relation in Eq. (73) is modified to 



which links in principle light and heavy neutrino parameters with triplet parame- 
ters. 

5.3. Left-right symmetric theories 

Left-right (LR) symmetric theories are a popular and appealing extension of the 
Standard Model, in which SU (2) l x SU{2)ji x U {1)b-l is the extended gauge group. 
Such a gauge symmetry can be arranged in breaking patterns of larger groups such 
as S'O(IO) or the Pati-Salam group. It is a natural framework to justify the type 
I + II seesaw terms in Eq. (20). The Higgs sector of the theory contains^^^ a 
"left-handed triplet" A^ with quantum numbers (3, 1,2), a "right-handed triplet" 




(76) 
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An with quantum numbers (1, 3, 2), and a bi-doublet $ with (2, 2, 0), which is not 
important for Quj3j3. At fow energies the potential consequences of LR symmetry are 
mainly right-handed currents mediated by a Wr with couphng strength gL = gn — 
g, and the presence of and A/j. Here A^ can be responsible for a contribution 
Ml = hvL to neutrino mass. It can give a direct contribution to Qu^jS, as discussed 
in Section 5.2, where we have learned that it is suppressed. The A^j gives mass to 
the right-handed neutrinos Mj,; = f v^, where vj^ is the vev of its neutral component 
and / a Yukawa coupling matrix. Often is is assumed that a discrete LR symmetry 
holds in addition, in which case h = /*, or AIl ~ M'^. Consequently', with writing 
mo ^ y V it follows 

T-l f '^''^ T -^ \ ( ^^T-l\ 
mv = ML-moM„ mo ^ vl { h y f y ] ^ vl { h y h y). 

\ VRVL J \ VRVL J 

(77) 

From the analysis of the scalar potential it follows vl oa jv-R and therefore neu- 
trino mass is zero in the limit ii^j — > oo, in which case there are no RH currents, 
because M-w^^ ~ gRVR. This connection of small neutrino mass and almost maximal 
parity violation makes LR symmetric theories very interesting. 

In what regards Oi^PfS, there are now several diagrams which allow for it. In 
certain variants of LR symmetric models one of the diagrams will dominate over 
the other, but we will not enter discussion of the details, and simply give the limits 
arising from each diagram individually. 

First of all, the A^ can mediate the process in analogy to the diagram in Fig. 25. 
It couples to the Wr instead of the W, and the two emitted electrons (as well as 
the quarks) are right-handed instead of left-handed. The amplitude goes as 



\MwaJ \MwnJ mj^^ 

where {AlR)ee can be written as with Mi the right-handed neutrino 

masses whose mass matrix Mr is diagonalizcd with V . Comparing with the naive 
amplitude in Eq. (68) gives ~ (m|^^ Af^^)/|(Affl)ee|, and from the standard 
amplitude (40) it follows 

|(M^)ee| < 10-15 GeV-^ (79) 



Expressing it with a dimensionless quantity is possible by defining 

|(Mii;)ee| '^Ip g 

'Often one considers /i = /, or = Mn, which happens when the discrete LR symmetry is 
connected to charge conjugation instead of parity. The Umits from LFV and CP violation in the 
quark sector case are stronger in this case'^^'^ . 
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Fig. 27. Left: quark level Feynman diagrams for left-right symmetric realizations of neutrino-less 
double beta decay with heavy neutrino exchange. Right: the corresponding diagram at LHC. 



The limit is compatible with TeV-scale left-right symmetry, one could for instance 
rewrite it as 

^/VsOOGcV^^/' 



Mw^ Z 1.9 



V 500 GeV / 



TeV. 



(81) 



In fact, limits from Oi^/3/3 arc competitive to other means of probing the parameters 
associated to LR symmctry^^^^^^^. Higgs triplets can be produced at the LHC, in 
the same way as shown in Fig. 26. Their decay into electrons or positrons probes 
fee = {Mji)ee/vn. An interesting possibility in these models is that Ml dominates 
the type I + II seesaw formula: rrii, = Ml- In this case Mj^ oc m^, i.e. the heavy 
neutrinos get diagonalized by the PMNS matrix and {M}i)ee becomes less arbitrary. 
However, in those cases it turns out that constraints from LFV, in particular ji — > 
3e, which can be mediated by triplets at tree level, force mA ^ Mi. This in turn 
implies that heavy neutrino exchange in connection with RH currents gives a larger 
contribution to Oj^/?/?^^*. 

Recall that heavy neutrino coupling to the usual LH currents is suppressed by 
small mixing m^i/Mij. The diagram to study is therefore the standard one from 
Fig. (8) with Wr exchange^^^, shown in Fig. 27. The amplitude goes as 



A 



Wn 



Elk 



(82) 



If / = /i (or f = h*) and type II dominance holds, then V ~ U {V ~ U*) and the 
PMNS matrix appears in this expression. By noting that the NMEs are the same 
as for the heavy neutrino exchange discussed in Section 5.1, we can use the limit 
from Eq. (71) to find 




r-5 



for ™Gc, 



for ^'Se, 



4.3) X 10-16 GcV 
16.6) X 10-16 GeV"^ 
16.6) X 10-16 GeV-^ for ^°°Mo 
4.3) X 10-16 GeV 



(83) 



fori^OTe, 
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Fig. 28. Quark level Feynman diagrams for left-right symmetric realizations of neutrino-less double 
beta decay. Left is the A-mechanism, right the r)-mechanism. 



and hence the dimensionless particle physics parameter has the same Umit as T^h in 
Eq. (72): 



mw 
Mwr 



< 1.7 X 10"* , 



(84) 



This hmit again corresponds to TeV scale, which can be seen by rewriting it as 

^ , , /500GeVV/' 
Mwn Z 1.5 I 1 



TeV. 



(85) 



A straightforward phenomenological LHC aspect of heavy neutrino exchange in 
left-right symmetric theories is seen in Fig. 27. Like-sign lepton production^^^ is 
possible and allows to directly test this mechanism. The current limit on Mwr set 
by LHC data is 1.4 TeV, both for very light right-handed neutrino mass^^'', as well 
as for masses between^^® 100 GeV and Mvi/„. In the future, LHC can detect masses 
up to a few TeV, and right-handed neutrinos up to TeV. This will test contributions 
of right-handed neutrino exchange to Qv^jS. 

The remaining two diagrams for LR symmetry are stemming from mixing of 
the left- and right-handed sectors. First of all, one of the W bosons in the standard 
diagram could be right-handed, leading to the left diagram in Fig. 28. Its amplitude 
is 



(86) 



where S is the matrix which quantifies the mixing of the SM leptons with RH 
currents. Note that one of the hadronic currents is right-handed. The dependence 
on 1/q can be understood from the RH nature of one of the vertices (see the 
comments after Eq. (36)). The dimensionless particle physics parameter is 



(A) ^ 77A = 



\Mwr 



(87) 
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The other contribution takes W-Wn mixing, quantified by tan^, into account and 
has an amphtude given by 

Ar,-Gltl^nCY,Ue^Se^^. (88) 

Here both hadronic currents are left-handed. The dimensionless particle physics 
parameter is 

(r/) ee; = tanC \ j2 -^^^ ■ (89) 

Note that the usual way in which both diagrams in Fig. 28 are drawn may be 
confusing. They are actually long-range diagrams with light neutrino exchange, 
and the lower vertex receives due to mixing with the RH current a (small) factor 
Sei — mo/Mji. This term requires non-zero and less than infinitely heavy. 
Without Mn and hence without lepton number violation it would obviously not 
be there. The implicit lepton number violation necessary for the existence of the 
diagram is illustrated by a Majorana mass term and a Dirac mass term, which gives 
a total contribution mu/Mn. 

In both the A and the rj diagrams one of the emitted electrons is right-handed. 
The nuclear physics becomes more complicated now, because the momentum de- 
pendence of the amplitudes, A oc = (uj,q), which introduces matrix elements 
corresponding to the time and space components of . In fact, the space com- 
ponents can give rise to 0^ — > 2+ transitions, whose observation would therefore 
be a clear signal*^ of the presence of right-handed currents in 01/(3/3. The main 
point here is that the time component (w) parts turn out to be suppressed by order 
{El — i?2)/w ^ 10^^ due to cancellation of the two diagrams with interchanged 
electron lines, whose energies are Ei and i?2. This suppression makes the time 
component parts of the same order as the space component {q} parts. The latter 
contain for (77) (not for (A)) two extra matrix elements, one of which stems from 
the nuclear recoil Q ^ q (with the electrons emitted as s-wavcs). This contribu- 
tion dominates and compensates the {Ei — £'2)/^ suppression. These features are 
explained in detail e.g. in Ref.'^. 

We are not aware of any recent comparative study of the relevant NMEs for 
these processes. Rcf.^° has recently summarized the calculation from^^^. We add to 
these results the ones from^^° (which do not contain ^^°Nd) and take the two sets 
of calculations as a span of NMEs. The result is 

(4.0- 11) X 10-^ from^^Ge, 

(1.4-4.4) X 10-^ from**^Se, 

(5.4 - 100.6) X 10^9 from ^°°Mo , , . 

(4.0-6.2) X 10-9 from"°Te, ^ ' 

(1.2 - 1.6) X 10"*^ from ^^^Xe , 



iv) = Vn < ' 



1.4x10-*^ from^^^Nd. 



The best limit from ^■^'^Tc of about 6x10 ^ corresponds to the naive result obtained 
by comparing the standard amplitude (40) with Arj, from which (?/) ^ 5 x 10"^ 
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is found. Again, comparing on the particle physics amphtude works amazingly 
well. However, this is here somewhat accidental, because the usually dominating 
time component part is suppressed by a factor 10~^. This suppression in turn is 
compensated by terms from the space components such as the nuclear recoil term. 

However, for (A) one would be two orders of magnitude off. While the naive 
estimate would give (A) ^ 5 x 10^^, the more correct procedure described above 
gives 



(A) = r/A < < 



(6.1 - 15) X 10"^ from^^Ge, 

(1.8- 3.8) X 10-*^ from^^Se, 

(9.8 - 54.5) X 10-^ from ^°°Mo , 

(5.8 - 8.9) X 10-^ from "°Te , ^ ' 

(2.1 - 3.2) X 10-6 from "^Xe , 

1.4x10-*^ from^^^Nd. 



Again ^-^^Te dominates the constraints and gives (A) ^ 9 x 10^^. The two orders of 
magnitude difference with respect to the naive limit originate from the suppression 
of the dominating time component part of the amplitude by a factor of electron 
energy ~ MeV divided by neutrino momentum ^100 MeV. 

To sum up. the full glory of left-right symmetric theories provides several possi- 
ble diagrams for Ovfi^: standard, heavy neutrino exchange, heavy neutrino exchange 
with RH currents, left-handed triplet, right-handed triplet, A and 77. In principle, all 
should be considered at the same time, yielding correlated constraints^^^'^^^'^^^'^^^ 
in a multi-dimensional parameter space spanned by parameters Mwr, tan^, m/^^, 
{Mfi)ee, J2 ^ei/^^i ^^'^ J2 UdSei- Ouc Can cxpcct that left-handed triplet and heavy 
neutrino exchange with LH currents can be neglected, but the remaining diagrams 
could give observable Oj^/3/3 if the relevant masses and scales do not exceed TeV 
too much. These scales correspond to values testable at the LHC, via lepton flavor 
violation or rare processes in the quark sector, and interesting works analyzing this 
interplay have recently been publishcd^^^'^^"*. 



5.4. Super symmetric theories 

In the context of supersymmetric theories i?-parity often plays an important role. 
It is defined as (_ i)3S-i-l-i-2s^ where B (L) is baryon (lepton) number and s spin. 
For particles R ~ 1 while for sparticlcs R — —1. The usual MSSM Lagrangian^^'^ 
conserves R. If R is violated, the following renormalizable and gauge invariant 
Lagrangian is allowed: 

Cfi = \jk U Lj el + A^^fc L, Q, 4 -f X^^^v^ + e, L,H^ . (92) 

Here the Li (Qi) are supcrficlds which contain the SM lepton (quark) doublets 
as well as the corresponding slepton (squark) doublets, uf, d^, e? are superfields 
containing the singlets of particles and sparticlcs, while contains the Higgs and 
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Higgsino doublets; i = 1,2,3 is the family index. The terms proportional to A"^^ 
violate baryon number, and have to be tiny to forbid too rapid proton decay. The 
remaining three terms in Eq. (92) violate lepton number by one unit. If i?-parity is 
broken™ (see^^® for a review), diagrams with two such vertices can therefore lead 
to neutrino- less double beta decay^^^^^^"'^"^^. Note that now there are two AL = 1 
vertices instead of one explicit AL = 2 mass term. 

There is again an interesting possible interplay here, namely that the i?-parity 
violating (RPV) terms are responsible for the neutrino mass itself. It is well-known 
that loop-induced Majorana neutrino masses can be generated by the A and A' 
terms. This would be their indirect contribution to Oi'/?/?, to be compared with 
the direct contributions to be discussed in this Section. A systematic analysis of 
the interplay of direct and indirect contributions to QvjSp is still lacking, but most 
often the neutrino mass constraints are weaker than the ones from Oi^/3/?, or the 
parameter space is chosen such that the RPV contributions to OvPfi dominate. 

The most simple possibility here is "bilinear i?-parity violation" , in which only 
the term LiHu is present. This is a realization of the type I seesaw, with the 
Higgsino playing the role of a single (TeV scale) heavy neutrino, which means that 
only one light neutrino is massive. Radiative corrections can generate the necessary 
other neutrino niasses^'^^. Bilinear i?-parity violation leads to mixing of neutrinos 
with neutralinos and of charged leptons with charginos. Effective cIl-u-c'^, e^-e'^- 
v, dfi-d-v and ul-u'^-v vertices arise, and can lead to Ov/S/S in diagrams^^"* which 
are similar to the ones in Fig. 29. For instance, one could have the upper left 
diagram with W instead of cl exchange, or the lower left one with neutrino and W 
exchange, instead of the and x, respectively. Those diagrams have been found to 
be suppressed with respect to the standard mass mcchanism^'^'^. We will concentrate 
on the trilinear terms from now on. 

The two RPV contributions to Qi'/SfH are shown in Figs. 29 (neutralino/gluino 
exchange, or X[-^^ mechanism) and 30 (squark exchange, or X[^^ X[^^ mechanism). 
Here the pion exchange dominance®'^'^® mentioned in Section 3.2 is realized. Limits 
on RPV SUSY parameters from Oi^/3/3 have been derived in^34,235^ 

The short distance diagrams are shown in Fig. 30. The naive estimate for the 
amplitude is 

- 1^ ' (93) 

^^SUSY 

where we set all sparticle masses to the same SUSY scale AgusY and the only 
relevant coupling is X'm, because the other vertices are order one gauge couplings. 
Comparing with the standard amplitude (40) gives Af^i/AgusY ^ '''x 10~^^ GcV~^. 

™In principle also the case of i?,-parity conservation can via box diagrams lead to 0^/3/3, if L 
violating sneutrino mass terms are present^^*. These are connected to the amplitude of Ou/B/S and 
L violating Majorana neutrino mass terms in analogy to the black-box theorem: the presence of 
one of the three implies the presence of the other two'^'^^'^^^ . However, the constraints from the 
sneutrino contribution to neutrino mass are stronger than the ones from 
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Fig. 29. Quark level Feynman diagrams for short-range i?-parity violating SUSY contributions to 
OuPP, which are proportional to 



Fig. 30. Quark level Feynman diagram for long-range /?-parity violating SUSY contribution to 
Ovfifi, which is proportional to \'^^^ -^113- 



The scale Agygy differs for the six diagrams in Fig. 29, for instance it is related to 

to|^ in the upper left and rrig in the upper right, etc. 

Note that in the diagrams x denotes all four neutralinos, which are linear com- 
binations of neutral gauginos and Higgsinos. In case that gluinos and/or squarks 
are exchanged, Fierz transformations have to be performed to obtain colorless op- 
erators. As a result'', scalar, pseudoscalar and tensor matrix elements arise. At the 
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end of the day, one can express the matrix element as^^^'*^'^'^^ 



(94) 



where Mf = 0(100), Mf = C(10) arc 2 nucleon NMEs and W = O(IOO) pion 
exchange NMEs (their absohitc magnitude exceeds the one of A^?^). Whether 
neutrahno or gluino exchange dominates depends on the SUSY parameters. 
The particle physics parameters in Eq. (94) are^'^"'^'^®'^^^ 



TTQfs X'lii nip 111 1 



'^'^ 2 Gj, ^ m., I m| + mi m? ) ' 

■fr i=l X' \ "I, dfi f-L dji / 

^ 27ra3 Afii 1 



Afii ^ TTip V^^(e) 
'^F ^ "^x. "^i 
Afn A mp /yL.(w)l/fl.(d) 14.(?/)1/L.(e) 14. (e)Ffl. (d) ~ 



/Viii v-^ nip 



where a^, oli are the S'J7(3)c and SU{2)l fine structure constants, respectively, and 
y are rotation matrices to go from the gaugino/Higgsino basis to the neutralino 
basis. As an example, consider gluino and pion exchange dominance, in which case 
the product of matrix elements and particle physics parameters in Eqs. (94,95) 
simplifies to 

-^#^fl+f^)T^^- (96) 
P'l 6 Gj, mgmj^ y \m.u^ J J 

The relevant NMEs are^^ between 387 - 569 for ^^Qe, 375 - 594 for ^^Sc, 412 - 589 
for ^°°Mo and 385 - 540 for ^^°Te. Hence, the current limits on 77^ are 

(4.9- 7.5) X 10-^ for^^Ge, 

(0.9-1.8)xl0-« for^^Se, 

- 1 (0.8-1.1) X 10-« for i"°Mo, 

(5.5-7.7) X 10-^ for "°Te. 

The value of £ 7.5 x 10~^ translates into Afii/(mg (1+mj /to|^)^) < 1.8x 

10""'^^ GeV~^, in very good agreement with the naive limit Aj^j^j/Agugy SS 7 x 10~^* 
GeV-5. 

Of course, supersymmetric particles are expected to be produced at the LHC, 
and Refs.^^^'^'^^ have recently analyzed the interplay of RPV contributions to 
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Fig. 31. Left: resonant seloctron production at the LHC as test of the short-range Ajj^j RPV 
diagrams in Fig. 29. Right: B'^-B^ mixing as test of the long-range A'j^gj^ ^'113 RPV diagram in 
Fig. 30. 




300 400 600 BOO lOQO 1 .9e+25 le+26 le+27 

observed Ge half life (yrs) 

Fig. 32. Left: mSUGRA parameter space (mo vs. ™ which single slcpton production may be 

observed at the LHC with y/s = 14 TeV and 10 fb~^ of integrated luminosity. The labelled contours 
show the search reach given by the labelled value of A^jj^. The white, dark-shaded and light-shaded 
regions show for '^^Ge that observation of single slepton production at the 5cr level would imply 
TIJ^ < 1.9 X 10^5 yrs, lO^'^ yrs > T^J^ > 1.9 X 10^^ yrs and T^J^ > 10^^ yrs, respectively. 
The upper and lower dashed curves show where the contour between the dark-shaded and light- 
shaded regions would move to if (mee) =0.05 eV were included with constructive or destructive 
interference, respectively. Right: ratio of the RPV amplitude (94) and the total amplitude of 
01//3/3 vs. the half-life if mo = 680 GeV and mi/2 = 440 GeV. Taken from^^S. 

0^(3/3 and collider physics. In particular, resonant selectron production^^^, ud'^ ^ 
CL — > ex — > euuL — >■ eued'^, was studied. The Feynman diagram is sketched in 
Fig. 31; note the typical like-sign dilepton structure. The first and last reactions in 
the chain involve X'm and the parton level cross section is proportional to Xfn/s. 
A numerical scan of a mSUGRA-like breaking scenario with mo and m.1/2 between 
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40 and 10"^ GeV, vanishing trihnear couphng Aq and tan/3 = 10 was performed, 
the result of which is given in Fig. 32. In the white region, resonant selectron pro- 
duction is forbidden by the current Oi//3/3- limits on ''^Ge. Observation in the darker 
shaded region implies that Qvfip should be observed in GERDA. Hence, if Oi//3/3 is 
discovered, searching for resonant selectron production at LHC is a direct test of the 
X'lii hypothesis. In light-shaded regions one does not expect observation of Oz^/3/3, 
and would hence rule out a possible contribution to the process. 

There is the possibility that the i?-parity violating diagram and the standard 
one contribute simultaneously (see Section 6). A possible effect of this is shown in 
Fig. 32: constructive interference of {mee) = 0.05 eV would move the interesting 
dark-shaded region up, and render observation of resonant selectron production 
very difficult. Destructive interference would move it down and make it easier. The 
right plot in Fig. 32 shows the ratio of the i?-parity violating amplitude Eq. (94) 
and the total amplitude of Oi^/3/3, for a particular point in parameter space. Ex- 
tracting the value of X'm from LHC and measuring the half-life of Qvfiji fixes this 
value. 



The long-range diagram from Fig. 30, given first in^^^^, involves no suppression 
by neutrino mass, and the amplitude can be estimated as 

-^k -Gf- Uei -3 Ksi Ki3 ■ (98) 

Here we have set all SUSY masses to a common scale Asusy, and took into account 
that the b-b'^ mixing is proportional to m;,/AsuSY (see below). Comparing with the 
standard amplitude Eq. (40) gives the constraint X'l^,! '^'113/^susY ^ lO^^"' GeV~'^. 
A more precise calculation gives constraints on the following quantity 

The angle 9^ and the masses to? in this expression arise from diagonahzation of 

01,2 

the symmetric matrix 

„.^^ -0.42 cos 2/3 -TO6(A + Aitan/3) \ 

. + ^2 _ 0.08 Af| cos 2/3 j ' ^^^> 

where tan /3 is the ratio of up- and down-type Higgs vevs, jj, is the /i-parameter, 
the trilinear coupling of Higgs scalars and fermions, and to? (to? ) the soft masses 
of the SUSY partners of the left-handed (right-handed) b quark. Nuclear physics 
is again dominated by pion exchange^'^^, with the relevant NMEs 2 to 3 orders of 
magnitude larger than the 2 nucleon NMEs. The spread of NMEs in^^ is 396 - 728 



\ 01 02 
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for ^'^Ge, 379 - 720 for ^^^Se, 405 - 691 for i"°Mo and 382 - 641 for ^^°Te. One finds 

f (4.0-7.3) X 10-^ for^^Ge, 



The agreement with the naive limit is very good. We have only considered here the 
b squark diagram. There are identical diagrams with d and s squark mixing, propor- 
tional to rrici X'lii X'lii and X'121 '^'1121 respectively. The first case depends there- 
fore on the same parameters as the neutrahno/gluino diagrams discussed above, 
but due to its dependence on rrid it is suppressed. The diagram with s squark mix- 
ing can be shown to be sub- leading due to strong limits from K^-K'^ mixing^'*'^, 
in which at tree level sneutrino exchange takes place. Those limits are of order 
\'-^2i ^'112 ~ 10~^ (Asusy/100 GeV)^, whose dependence on the parameters is easy 
to understand. About the same order are the limits on X'n^ from B^-B^ mixing 
(Fig. 31 sketches the relevant Feynman diagram), which have to be compared with 
X'l^i -^'113 ^ 10~^ (Asusy/100 GcV)'^ from neutrino-less double beta decay. This im- 
plies an interesting interplay of B physics and Oi/ (3/3: as long as the SUSY breaking 
scale docs not exceed TcV, the limits are similar. However, as the B'^-B'^ mixing 
diagram proceeds with sneutrino exchange and the Oi'/J/S-diagram with b squarks, 
a more detailed analysis is in order, which has been performed in Ref.^^^. As a 
result, the constraint is currently stronger than the one from Oi^/3/3, but can 

be responsible for observable half-lifes of 10^^ - 10^^ yrs for '^^Ge, which was the 
isotope studied in^^^^. In analogy to the right plot of Fig. 32 one could again define 
a ratio of matrix elements and study its range as a function of the half-life^'^^ . 

5.5. Majorons 

The term Majoron denotes very light or massless particles x*^ which can couple 
to neutrinos. Originally Majorons were Goldstone bosons of spontaneously broken 
global lepton number. This Majoron could be part of a weak singlet^'*^, doublet or 
triplet'^^^ , the latter two cases being ruled out by their unacceptable contribution to 
the Z width. Another set of important constraints stems from astrophysics^*'^'^''^. 
In the context of triplet Majorons it has been noted that the decay mode^''^^ 



is induced; see Fig. 33. Several different approaches of (almost) massless scalar 
particles coupling to neutrinos and their impact on Oj^/3/3 have been made in the 



past246'247,248,249,250,25i,252,253,254^ rpj^^^^^ include sccnarios in which Majorons are 



not Goldstone bosons, or carry lepton number, such that lepton number is actually 
conserved and Qi^/3(3 is forbidden. Other examples are when Majorons are vector 
particles^^"'^, or doublet Majorons^^® in which the Majoron is the SUSY partner of 
the neutrino. Extra-dimensional Majorons with a set of Kaluza-Klcin modes was 




(101) 



iA,Z) {A,Z + 2) + 2e- +x' 



,0 



(102) 
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Fig. 33. Quark level Fcynman diagram for a one and two Majoron realizations of netitrino-Iess 
double beta decay. 



also proposed^^^. Simple singlet Majoron models allow coupling of x° to right- 
handed neutrinos with strength m^/M, where M is the scale of spontaneous lepton 
number breaking, hence M ~ M^. Thus one does not expect sizable coupling and 
0://3/3-rates. This is different in more complicated models. It was also realized that 
decays with two Majorons in the final state are possible^^^: 

(^,Z)^(A,Z + 2) + 2e- + 2x°, (103) 

see Fig. 33. At the end of the day, one writes the decay rate as 

T'''' = \{g^)f°''''\M^\'G^{Q,Z), (104) 

where \ {gx) \ is an averaged and model-dependent coupling constant, its power ob- 
viously depending on single or double emission. The phase space factor Gx{Q, Z) 
depends also on the number of final state particles, but also on the model, in par- 
ticular on the nature of the Majoron. The experimental quantity to distinguish 
Majoron modes from Oj^/3/3 is of course the energy spectrum of the two emitted 
electrons. In the original triplet model, with Qx^x^ the coupling of the Majoron 
with two neutrinos, the amplitude can be written as ^ ~ G^p (ffx)/?^' which has 
one dimension of energy less than the amplitudes considered before, because the 
phase space integration for one additional final state particle implies two powers of 
energy. Hence, the decay width goes as Q*" {n = 1) instead of for QvjSfi. Similar 
models with double Majoron emission have consequently a decay width propor- 
tional to [n = 3). With the same logic it follows that single Majoron decays 
where the coupling goes with d^x have a width proportional to {n = 3), while 
double Majoron decays go with Q^^ {n = 7). The integer number n in the above 
considerations indicates the "spectral index" of the two electron spectrum^^^ 



dEi dE2 « - ^1 - ^2)" - < V ^2 - ml El E2 , (105) 

neglecting Fermi functions and prefactors. Reasonable estimates could now be 
made, again by comparing the amplitudes, and taking into account the different 
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Table 13. Categories of Majoron models as first proposed in^^^. 
Given are whether OvPfi can take place, if single or double Majoron 
emission is predicted, the spectral index and the lepton number 
of the Majoron, whether it is a Goldstone boson, and the limit on 
its coupling, taken from'^^'^. The limit for IF is estimated here. 
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1.7 X lO-"* 
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/ 
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IE 


/ 
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1.5 
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x° 


2 
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~ lO-'' * 


IIB 






1 


-2 




1.7 X lO"'' 


lie 




x° 
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-2 
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0.024 


IID 




x°x° 
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-1 




1.5 


HE 




x°x° 


7 


-1 


/ 


1.3 


IIF (vector) 




x° 


3 


-2 




0.024 



Note: *this is a limit on g^l^ jM in units of GcV"'^, where g is 
the ■)^vv coupling and M the low energy string scale in the extra- 
dimensional framework studied in^^*. 



phase space dependence and a factor 2(27r)'^ for each additional phase space inte- 
gration. In this way one finds for instance that for single Majoron modes with 
n = \ the standard contribution [0"]^ {rriee) / (f'Y' has to be compared with 
(Gf {9x)/l^y <97(2(27r)3), from which it follows (g^) < 10-^ and forn = 3 that 
((/^) ^ 1. Nuclear physics aspects are dealt with in255,256^ g^^^^ ^j^l not go into 
detail here. For single Majoron and n = 1 cases the NMEs fr-om the standard 
interpretation can be used, while for the other cases different NMEs need to be 
calculated, similar to the situation for the (A) and (rj) terms in the presence of 
right-handed currents, discussed in Section 5.3. We rather summarize the limits on 
the various model categories, which first have been described in^^^. This is shown 
in Table 13. 

5.6. Other mechanisms 

We will discuss other proposed realizations of Oiy/3/3 in this Section. 

Non-renormalizable effective operators in the Lagrangian Ccs = 

Q4:+d^jt^d ^g^j^ generate neutrino Majorana masses and/or lepton number violation, 
the most simple example being the Weinberg operator of dimension 4 + d = A + 

I = 5 in Eq. (17). Operators with AL = 2 have been classified up to dimension 

II in258,259^ They can generate neutrino mass directly (the Weinberg operator) 
or via loops, by closing some of the external legs. It is also possible that those 
lepton number violating operators generate a direct contribution to OvfS/S (note 
that is effectively a ududee operator, which has dimension 9). There are 
five dimension 9 and fifteen dimension 11 operators which have this property^^*. 
One example is = LLQQd'^d'^ . Closing the external Q and d'^ lines with Higgs 
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loops gives a neutrino mass term of order rrii, ~ m^/A/(167r^)^. A limit on A is 
estimated from the direct contribution of to which has an amplitude of 

order A^^. Hence one limits A <; 3 TeV, and therefore a tiny mass of rrii, ^ 10~^ 
eV is generated. All dimension 9 operators generate a limit of order 3 TeV on 
their associated suppression scale. Regarding dimension 11 operators, their 0^/3/3- 
amplitude can be estimated as w^/A^, hence A TeV. In^'^^ all 129 operators up 
to dimension 11 have been studied and the scale A has been fixed by requiring the 
operator to generate ~ 0.05 eV. This fixes their contribution to Oz^/3/3. Some of 
the operators can now be disfavored, because their direct contribution to Oz^/3/3 can 
be too large^^^. 

Leptoquarks can couple to quarks and leptons and the SM Higgs doublet, and 
have the potential to lead to lepton number violation and Oz^/3/3. Their properties 
are similar to i?-parity violating mechanisms of Qvfip. In Ref.^^" the vertices for 
{S, V^yd-v, (5*, y )-d-e, (5, and {S, V'^)-u-e interactions have been worked 

out, where 5* (F^) are scalar (vector) leptoquarks with electric charge — i or |. 
Effective u-e-v-d vertices arise and the coefficients depend on the original leptoquark 
couplings and masses, the latter obviously as M^y. Writing the amplitude naively 
as vAlq ~ Gp a/Mg y/q, where a the typical coefficient for the effective vertex, 
one finds limits oi a ^ 10~^ for 100 GcV leptoquarks, which is within one order 
of magnitude to the actual limits derived in^^°. In that paper the definition of the 
coefficients in terms of original parameters can be found. 

In^^^ scalar bilinears (coupling to two fermions) have been considered, and 
typically one dimcnsionful coupling /x and 3 propagators are present in the 0vf3/3- 
diagrams, leading to amplitudes of the form /i/Af^, where M is the common mass 
of the bilinears (see also^^^). 

Rather surprisingly, given the popularity of scenarios with extra spatial di- 
mensions, there are only few papers discussing its consequences on Qvfip. Ref.^^'^ 
showed that within ADD scenarios small Majorana neutrino masses can result if 
lepton number is broken on distant branes, with the breaking being communicated 
to our brane by messenger particles. Ref.^^"* used this finding and translated limits 
on the Majorana mass {m.ee) in limits on the number of extra dimensions, compact- 
ification radius of the extra dimension and messenger mass. A generic feature of 
extra dimensional theories is the presence of Kaluza-Klein (KK) excitations of par- 
ticles which feel the extra dimensions. If Majorana neutrinos do so, the associated 
tower contributes in principle to QvjSp. The case of all excitations being Majo- 
rana neutrinos was discussed within a particular model in^^^. Excitations heavier 
than 100 MeV will have NMEs with the characteristic features of heavy neutrino 
exchange discussed in Section 5.1. In the model considered in^^'^ two parameters 
had to be chosen, the radius R of the extra dimension and the brain shift parame- 
ter a, introduced to make the neutrinos with opposite CP parity couple to the W 
bosons with unequal strength. Constraints on those parameters are possible. Ref.^^^ 
studied an extra-dimensional scenario based on a warped Randall-Sundrum model 
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leading to low scale seesaw, in which also a tower of order GcV sterile neutrinos 
is present, which can be constrained by Models in which only gauge bosons 

or Higgs scalars possess KK excitations, such as in^^^, are dominated by the usual 
light neutrino mass mechanism. 

Scalar octet seesaw has been proposed in Ref.^^* to have TeV scale neutrino 
mass generation with sizable LHC cross sections. In this mechanism a one-loop 
diagram including a weak scalar triplet S and weak fermion singlets or triplets 
Pi, all octets under SU{3)c, produces a Majorana neutrino mass. The p and S 
particles could mediate double beta decay via the usual standard diagram with the 
W replaced by S and the neutrinos replaced by pi. The amplitude is proportional to 
c,^^/(TOg rup.), where Cud is the coupling of to u and d quarks, which is constrained 
from flavor violating transitions. 

A fourth generation Majorana neutrino with mass M4 behaves exactly as 
a heavy neutrino discussed in Section 5.1. Thereforc'^^^, it receives the constraint 
15641^/^1^4 < 1.8x 10~® GeV~^, with being its mixing with the electron. Pushing 
its mass down to collider level would require cancellation with other contributions 
to OujSp. 

Composite neutrinos^^° lead to heavy neutrinos A^* which are excited states 
corresponding to a scale of the SM neutrinos. Their coupling with gauge bosons 
goes with //Ac, / being a coupling constant, and the amplitude for Ovfi/S goes with 
P / Kc/Mm* and is sensitive to TeV scale exciteness^ ""^ . 

So-called 3-3-1 models with an initial S'L/(3)l gauge symmetry contain new 
gauge bosons and scalars, which can contribute to Oiy/3/3. Those cases have been 
studied in Refs.^^^^^^'^'^^'^, and constraints on the masses and mixings with the SM 
fermions have been obtained. Majoron emission is also possible in those models, 
because typically neutral scalars with lepton number exists, whose vevs induce 
spontaneous violation of lepton number, see Section 5.5. 

We conclude this section with more exotic proposals. Effects of scalar unparticles 
in have been discussed in^^*, and an unusual model with colored scalars 

coupling to leptons and quarks, which can mediate Ovf3/3, in^''^. Recently it was 
proposed that a huge number of copies of SM particles exists^^^, which could solve 
the hierarchy problem and, if a permutation symmetry is added, explain also small 
neutrino masses. It was shown^^*" that this leads to basically vanishing amplitudes 
for Ov/Sp. 

6. Distinguishing mechanisms for neutrino-less double beta decay 

We have seen in the last two Sections that there are several well motivated frame- 
works in which observable neutrino-less double beta decay can be expected. Obvi- 
ously, means to distinguish the various possibilities are necessary. This is a common 
problem for all experiments looking for new physics, for instance lepton flavor viola- 
tion, where observation of, say, /i — > £7 alone does not prove the presence of super- 
symmetry, but could mean a lot of different things (Higgs triplets, extra dimensions, 
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non-unitary PMNS matrix, etc.). In this Section wc will classify three possible tests 
of the underlying mechanism of neutrino-less double beta decay. Mostly the domi- 
nance of one mechanism is assumed, but we will discuss the simultaneous presence 
of more than one mechanism as well. 

6.1. Distinguishing via effects in other observables 

This obvious possibility has been discussed at several occasions in the last Section. 
In particular within i?-parity violating SUSY and left-right symmetric theories TeV 
scale particles can lead to observable Oi//3/3. Production of such particles at the LHC 
is then a check of these mechanisms, in particular if the like-sign dilepton signature 
can be used, such as for heavy right-handed neutrino production^^^, Higgs triplet 
decays^°^, or resonant selectron production^^''. However, checks are also possible 
in processes in which lepton number is not violated, but instead quark or lepton 
flavor is not conserved^^^. To perform such studies, one can express the relevant 
processes in terms of effective operators suppressed by some high energy scale. The 
scales of flavor violation and lepton number violation could be different, but are 
related or even identical in some cases. The flavor parameters on which Oz//3/3 and 
flavor violating processes depend can also be different. 

Examples mentioned above are B^-B^ mixing induced by ^'n^ couplings^"^^. 
Note that here the parameters corresponding to flavor {X'l^i ^iiz) ^-re the same for 
Oi//3/3 and B^-B^ mixing, but different particles are involved: squarks in Qvfip and 
sneutrinos in B^-B'^ mixing. In left-right symmetric theories an important con- 
tribution to lepton flavor violation stems from Higgs triplet exchange, which can 
mediate — > 3e at tree level. Here the flavor physics parameters (also the ones for 
fjL — s- 67) are not directly related to the ones which govern Qvpp. If TeV scale physics 
generates Oi//3/3, and if no special flavor structures and not too different flavor and 
lepton number violating scales are present, one expects^^^ a ratio R ^ 10~^ of the 
rates for ^-e conversion in nuclei and ^ ej. Therefore, if _R ~ 10^^ is observed, 
the standard interpretation of light neutrino exchange in QvjSp is favored'^. 

6.2. Distinguishing via decay products 

We have seen that the Lorentz structure of the different mechanisms of Qvfiji can 
be different. This implies that energy and angular correlations of the two emitted 
electrons may be different'^'^^'*^. In particular the SuperNEMO experiment will be 
able to perform such measurements, because the set up of foils with Oz^/3/3-isotopes 
in a magnetic field allows tracking of the individual electrons, instead of "only" 
measuring their total energy. The design of the detector allows direct detection 
of two electrons from double beta decay by a tracking chamber and a calorimeter 
measuring individual energies and times-of- flight. In Ref.^^^ the collaboration has 

"Massive neutrinos imply lepton flavor violation in decays like — > £7 at an unobservably small 
level. 
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Fig. 34. Left: constraints at la on the model parameters from an observation of Ovf}fi of *^Se at 
half-life 10^''' yrs (outer blue elliptical area) and 10^® yrs (inner blue elliptical area). Adding the 
reconstruction of the angular (outer, lighter green) and energy difference (inner, darker green) 
distribution drastically shrinks the allowed parameter space. Right: adding information from the 
decay of ^^''Nd. In this example, 30% admixture of the A-mechanism is assumed. Taken from^'^^. 

simulated the potential discrimination power between the standard mechanism and 
the A-mechanism within left-right symmetry. The differential decay width can be 
written as^''^'''^° 

dV ( {1 — /3i f32 cos 6) standard mechanism ^incN 

dEi dE2dcose °^\{Ei- E2f (1 + Pi P2 cos 6) A mechanism ' ^ ' 

where £'1,2 are the kinetic energies of the electrons, /3i^2 their velocities and 9 
the angle between them. One can define an asymmetry Ag = — ) / ( + 
iV_), where N+ (iV_) is the number of events with Q > -n j2 (6 < tt/2). Another 
asymmetry \s Ae ~ {Ny — N^) / {Ny + N^) , where Ny (N^) is the number of events 
with El— E2 < Q/2 (El — E2 > (3/2), where Q is the energy release. Fig. 34 shows 
a result from^''^, where a 30% error on the NMEs and the simultaneous presence of 
the standard term and 30% admixture of the A-mechanism has been assumed. The 
energy difference distribution turns out to have a stronger discrimination power. 

Another aspect of identifying the Oi^/3/3-mechanism with the decay product is 
when the Majorons as additional particles are emitted, in which case the energy 
spectrum of the electrons is different from the 0vP(3- or the 2://3/3-spectrum, as 
discussed in Section 5.5. 

6.3. Distinguishing via nuclear physics 

We have not spent much attention on the nuclear physics details of the Oi^/3/3- 
mechanisms, and argued mainly on the particle physics amplitude level. However, 
there is nuclear physics involved, and if its uncertainties can be kept under con- 
trol, it could in fact be helpful®^'^°'^^ to disentangle the various mechanisms of 
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Fig. 35. Predictions for the inverse half-life for different NMEs, Oi^/S/S-mechanisms and isotopes. 
Taken from^*^. 



OvP/3. The use of multi-isotope determination of 0^/3/3 to test NME models for the 
standard mechanism was discussed in^"^-^^. However, it should be clear that this 
is a quite challenging task, and in particular requires that the spread of the NME 
calculations is not much larger than the experimental error on the half-life. Nev- 
ertheless, the same strategy can be applied to disentangle different mechanisms of 
Oi//3/3. Fig. 35 shows the result of Ref.^°, where different isotopes, NMEs and mech- 
anisms of Ov/Sp were compared. Those were the A and 77 diagrams within left-right 
symmetry, heavy neutrino exchange and i?-parity violating SUSY. The NMEs were 
two sets of QRPA calculations (with their parameters fitted to reproduce single 
beta decay and 2i'/3/3, respectively) and a shell model evaluation. A 10% error on 
the calculations was assumed. The individual parameters of lepton number viola- 
tion were chosen such that for ''^Ge the half-life is the same for all mechanisms. 
As can be seen from Fig. 35, for different isotopes there can be a significant spread 
of the half-lifes. By simulating sets of Oj^/3/3-rates it was estimated that 3 positive 
experimental results are required to pin down the mechanism of Quj3p, if a total 
(theoretical, systematical and statistical) uncertainty of 20% or less can be achieved. 
For 40% uncertainty four results would be necessary. Analyses in similar spirit can 
be found in Refs.^^'^^. Obviously, multi-isotope determination is here crucial. 

Another possibility to distinguish the mechanisms is the rate of the ground- 
state-to-ground-state transitions to the rates of decays into excited states®^'*^'®^'^^. 
The latter could be 0^ or 2^ . The decay into 0^ is experimentally easier to identify 
because two photons associated with the transition first to 2^ and then to the 
ground state are emitted. Transitions to 2^ states have higher sensitivity to right- 
handed currents^^, and observation with large rates would signal the dominance 
of these mechanisms, in particular the A-contribution. The experimental situation 
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is summarized in^'^*^. it is worth noting that 2i^/3f3 to Of excited states has been 
observed. Observation of Oiy/3/3 into ground states and excited states in the same 
experiment could be used to rule out the possibihty that an unidentified background 
peak mimics the Oi^/J/S-signal^*^ . The development of high granularity detectors, 
large enough Oi^/3/3-rates and precise nuclear physics is necessary to realize this 
consistency test. In the standard interpretation the decay to excited states occurs 
with a rate suppressed by a factor Rex — 10^ with respect to the 0+^ — > 0+^ 
transition. This factor is a combination of kinematics {{Q — i?ex)/Q)^> Q ~ -E'cx 
being the energy release to the excited state, and nuclear physics \Mg.s./Mcx\'^, 
and is sensitive to the mechanism of Qvfip. For instance, Ref.^^ has found for ^^Ge 
suppression factors of i?ox = 96,48 and 120 for the standard mechanism, heavy 
neutrino exchange and gluino exchange in i?-parity violating SUSY, respectively. 
For ^"^^Xc the factors arc 17, 38 and 153, while for ^™Mo the result was 17, 17 and 
59. These differences may be used to distinguish the mechanism, if one assumes 
the nuclear physics uncertainties to be under control. In this respect we compare 
the NMEs for ™Ge from Ref.^s (QRPA) with the ones from^" (IBM)°. For ™Ge 
the QRPA NMEs for ground state and Of transitions are 2.80 and 0.994, leading 
to a factor 7.93 in the relative half-lifes. The IBM NMEs arc 5.465 and 2.479, 
hence a factor 4.84. For i°°Mo the QRPA NMEs are 3.21 and 1.76, thus a ratio 
3.33. IBM gives NME values of 3.732 and 0.419, thus a ratio 21.26. Therefore, the 
notorious NME uncertainty will again be a problem of the procedure described here. 
Nevertheless, important information to the field would be added by observation of 
Qvf3j3 into excited states. 

Another possibility^^'^^^ to disentangle the mechanisms is the ratio of Qvfifi to 
0J//3+/3+, 0j/;3+EC or O^ECEC, see Eqs. (7,8,9). For instance^^, the ratio of the 
rates of 0i//3/3 of ^^Qe and Oi//3+/3+ of ^^^Cd are about 2087, 30435 and 1826 for 
the standard, the A- and the ?/-nicchanism. respectively. For Oz^/3+EC of ^"^Cd the 
ratios are 148, 12 and 217. The same comments on nuclear physics uncertainties as 
for excited states apply here, in addition to the problem of even lower rates. Recall 
however the possibility of resonant enhancement"'^ of Oi^ECEC. Double electron 
capture to excited states has recently been discussed as another way to distinguish 
mechanisms^^^. 

6.4. Simultaneous presence of several mechanisms 

We will now discuss aspects of simultaneous presence of more than one Qv^jS- 
mechanism. The different mechanisms can add coherently in the amplitude, see 
Eq. (66), and interference effects arc possible. However, at leading order only terms 
in which the hclicities of the emitted electrons are identical can interfere. The fact 
that helicity is not exactly equal to chirality for the emitted electrons with energy Ef. 

° Given the progress made in recent years it may not be fair to compare a 10 year old calculation 
with a very recent one. However, we are not aware of any recent QRPA re-evaluation and expect 
the ratios to be more stable than the NMEs themselves. 
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leads at the end to a phase space factor of the interference term suppressed^^'^^ by 
one order of magnitude (corresponding very roughly to [Ep/Q)"^), and interference 
is almost negligible. For instance, in the standard mechanism both electrons are left- 
handed, while in the A-mechanism one is right-handed, thus these two diagrams do 
not interfere. The chirality of the emitted electrons is indicated in the respective 
quark level Fcynman diagrams which are shown in this review. It is conceivable that 
destructive interference of several mechanisms leads to a vanishing rate of in 
one or more isotopes. Nuclear physics differences of the relevant NMEs could, but 
do not have to, lead to a non- vanishing rate in other isotopes^^. 

The procedure to deal with the presence of several mechanisms has been outlined 
jj^283,93^ Consider first the presence of two essentially non- interfering mechanisms, 
e.g. light and heavy right-handed neutrino exchange. If two experiments using dif- 
ferent isotopes have found evidence for Oi'/SfS, one has (see Eq. (4)) 

where the superscript a, b denotes the two isotopes and the subscripts I and Nji 
denote standard light neutrino exchange and heavy right-handed neutrino exchange 
with Wr instead of W in left-right symmetric theories, see Eq. (82). Solving for the 
particle physics parameters gives 



2 _ 

~ I Ay(a|2 I A/ff" |2 _ I A^&|2 I A^a |2 ' 

(108) 



2 _ \ Mt?/{Tlf^G'^)~\M\?/{T-f ^ G-) 



ma 



RecaU the present limits of |77i| ^ 9.8 x 10^*" and \riNR\ ^ 1-7 x 10^^ from 
Eqs. (67,84). Fig. 36, taken from^"^ shows an example solution of Eq. (107). Knowing 
the half-life of one isotope constrains the half-lifes of the other ones. 

Consider now two interfering diagrams, for instance the standard mechanism 
and gluino exchange within i?-parity violating SUSY, where the emitted electrons 
are both left-handed. There is an unknown phase between the two contributions, 
and the total half-life can be written as 



iT^/2)-' = G [\Mi\' h\' + \M%^ P P + 2\Mi\ \M%^ \ \m\ Iv^^ I cos ^) . (109) 

Obviously, three positive observations of 0^/3/3 in three different isotopes are re- 
quired in order to extract the three independent parameters |r?i|, \r]^^ \ and cos0. 

An example from^"^ is presented in Fig. 36. The current limit is jry^ | < 7.5 x 10^^, 
see Eq. (97). 

7. Alternative Processes to Neutrino-less Double Beta Decay 

The last section of this review deals shortly with alternative processes to Oi^/3(3, 
i.e. alternative probes of lepton number violation. The presence of AL = 2 can 
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Fig. 36. Left: solving the ease of non-interfering proeesses in Eq. (107) for the indieated hypo- 
thetical half-Iifes of ^"^Ge and ^^°Te. The blaek (red) lines are |r;ip (|'7iVfjP)i the spread be- 
tween the solid and dashed lines arises from varying the NMEs between 5.44 — 5.82 (4.18 — 4.70) 
and 411.5 - 264.9 (384.5 - 239.7) for '^'^Ge (""Te). The horizontal dashed hne corresponds to 
(mee) < 0.2 eV. The physical (positive) solutions for Iryip and Iry^VjjP are constrained within the 
solid (dashed) lines. Right: two interfering processes with constructive interference cos0 = 1 and 
fixed NMEs. The solid line is the dashed line is |r;^ |, the dash-dotted line is 2 | cos<j). 

The blue areas are forbidden. Taken from^'^. 



manifest itself in going from Icpton number L = to L = ±2, typical for a decay 
process, or from L = ±1 to L = =f1, typical for conversion processes. It could also 
be that \L\ — 1 goes to \L\ = 3, e.g in lepton decays or collisions with initial lep- 
tons. Finally, in like-sign lepton collisions or Oi^ECEC one could go from L = —2 
to i = 0. 

Neutrino oscillation probabilities arc not sensitive to the Majorana nature of 
neutrinos. However, in principle Va transitions are possible, whose probabili- 

ties are unfortunately suppressed by the factor (mi/E)'^, in analogy to the standard 
mechanism of Ov/Sf]. There are in principle differences between Dirac and Majorana 
neutrinos, for instance it is easy to show that Majorana neutrinos do not have a 
vector current. Again, in amplitudes the difference of Dirac and Majorana neutrinos 
due to the absence of vector currents for the latter goes with rrii/E. This annoying 
property is known as the^®* "practical Dirac-Majorana confusion theorem" . 

A recent review on the electromagnetic properties of neutrinos can be found 
jj^28.5 short, Majorana neutrinos cannot possess diagonal magnetic moments, 
i.e. {ve) — ^ i^e transitions would only be possible for Dirac neutrinos. This can be 
seen by looking at the magnetic moment operators^®^ Mq/j cr^i^ {vR)pF^^'^ for 
Dirac and ^lapT^a '^tiy {^l)^ -^^^ Majorana neutrinos. In e scattering exper- 
iments the helicity and flavor of the final state neutrino cannot be measured and 
there is no way to distinguish Dirac from Majorana in this way. One possibility 
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would be via spin flavor transitions in supcrnovae, in which a magnetic field trig- 
gers {i'e)L ~^ ii^n)R, with subsequent oscillation of the active {v^)r into {i'e)R- The 
usual Ve neutronization burst can be heavily affected by this effect^^^. While the 
SM extended with massive neutrinos does generate too small magnetic moments, 
H (X m^, in some extensions of the SM it would be possible to generate the required 
large magnetic moments^^^. 

We have mentioned already high energy tests of lepton number violation such as 
the like-sign dilepton signature of heavy Majorana neutrino production^^^^^^^, Higgs 
triplet decays^"^, or resonant selectron production^^^. One may wonder whether 
there are low energy processes, in analogy to Oi//3/3, which can probe the effective 
Majorana mass, maybe even without any nuclear physics complications. However, 
the {mi/EY suppression of the rate together with Avogadro's number Na render 
QvPP the only realistic probe. With order kg of a 0;//3/3-isotope one has order Na 
atoms, which compensates the Dirac/Majorana factor (irii/qY . In principle, there 
are decays like — > tt" e+e"*", which depend on the effective mass in the same 
way as Qv[3p does, and do not suffer from NME uncertainties. However, calculating 
the branching ratio yields^^^ 



to be compared with the experimental upper limit^^'^'^'^'^ of BR(i4r+ tt^ e+e+) < 
6.4 X 10^^'^. If it was possible to increase the number of charged kaons by 20 or- 
ders of magnitude, one could go for decays like^^^ K'^ — > tt^ ("neutrino- 
less double muon decay") and test (m^^), i.e. the other entries of the mass 
matrix^^^'^^'^'^^^'^^^'^^^'^^^. Other decays which have been studied in the past 
include lepton number violating decays of t leptons^^^'^^^. top quarks and W 
bosons^^^, D and B meson decays^^^'^^^''^™'"^*^^, or hyperons"^^^ . Collider processes 
such as^^^ v^N ^ XfjL~ P'^ or^°^ e~p — >■ XVf. (3'^ have also been discussed. In 
addition, searches for conversion processes such 

g^g3,305,306 ^ 

' {A,Z) e+{A,Z-2) 

QJ.307 ^- j-^^ _^ /^^(^, Z — 2) have been proposed. 

For very light and very heavy Majorana neutrinos the above processes are not 
very helpful. Recall however the general property of Majorana neutrino exchange 
as displayed in Fig. 23: for neutrinos whose masses correspond to the typical energy 
scale of the process the sensitivity is largest. Taking if + — !• tt" ijl^ as an example, 
Ref.^^^ has obtained very strong limits on masses between 245 MeV and 389 MeV, 
with |C/^ip down to the 10~^ regime. The constraints are strong because there can 
be "s-channcl" diagrams. Other decays have been analyzed in-^oo,308^ 

Not many works exist which study the above processes in non-standard mech- 
anisms. Examples include Ref.'^°^, where meson decays such as K'^ — 7r~ 
mediated by i?-parity violating SUSY were found to provide no significant limits. 
The same was shown in'^^" for (/i~,^"*") conversion, or for (/i~,e+) conversion in 
various mechanisms'^^ ^. Doubly charged Higgs exchange in — > 7r~ /i"*"/!"*" was 
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N 



(a) 



w- 





w- 




w- 



(c) 



w- 



Fig. 37. Fcynman diagram for inverse neutrino-less double beta deeay. Diagrams (a) and (b) are 
Majorana neutrino A'^ exchange, diagram (c) triplet exehange. 



also found to generate negligible rates 



312 



A particularly clean probe of lepton number violation is "inverse neutrino-less 
double beta decay" . This is not (A, Z + 2)++ + 2e~ ^ {A, Z), but 



(111) 



This reaction can be tested if a future linear collider is run in a basically background- 
free like-sign mode, and has frequently been proposed as a probe of LNV and new 
physics in gcncral^^''. The process does not involve any nuclear, hadronic or atomic 
uncertainties or difficulties and is presumably the cleanest probe of lepton number 
violation. If Majorana neutrinos are exchanged, sec Fig. 37, and with neglecting 
the mass of the W , the cross section reads 



da 



dco& 9 



t 



(112) 



where t and u arc the usual Mandclstam variables, {M^)i is the mass of the neutrinos 
(including light and heavy Mi) and Vei their mixing with electrons {N^i and 
Sei)- There are interesting special cases for the cross section: 

• if only light active Majorana neutrinos contribute to the process, then the cross 
section is 



cr(e e 



< 4.2 X 10" 



Gp , ,2 
4 7^2 



leV 



(113) 



fb. 
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Fig. 38. Cross section for e'e' W'W- with = 1 TeV (left) and = 4 TeV (right) and 
three Hmits for the mixing parameter IVeip. The dotted line corresponds to five events for an 
assumed luminosity of 80 (s/TeV^) fb~^. 



hence far too small to be observable; 

• if only heavy Majorana neutrinos contribute to the process, then we can bound 
the cross section using the Oj^/3/3-limit from Eq. (71) as 

(114) 

again far too small to be observable; 

• the high energy limit of -^/s — >■ c» is 

a{e-e- ^ W-W') = ^ (Af.).)' ■ (115) 

This seems to violate unitarity, because the cross section for an s-wave process 
should vanish in the high energy limit. However, recall the exact seesaw relation 
^7Vg-mi + Sl^Mi = 0, as discussed in in Eq. (73). This relation guarantees 
that the cross section vanishes in the high energy limit. 

While small and large masses cannot give sizable cross sections, intermediate scale 
neutrino masses {Mi,)i ~ -^5 can give appreciable event numbers, as expected from 
the general behavior of LNV processes with Majorana neutrinos involved. In Fig. 38 
we show an example of the cross section as a function of neutrino mass. Different 
limits on the mixing Vei are inserted: no limit, the global limit and the limit as 
implied from the limit on (^) = \Vf,i\^ /M^. Note how the case of iKiP = 1 follows 
the general trend of Fig. 23. The above processes can also be searched for at e~ ii~ 
or fJ-^ fi^ machines. 

The process e~e~ W~W~ can also be mediated by a Higgs triplet, but due 
to small neutrino masses has a tiny cross sections unless a very narrow resonance 
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is met. One can show that if neutrino and triplet exchange occur simultaneously, 
the unitarity of the cross section is saved by the relation in Eq. (76). Right-handed 
Higgs triplets or Wr could be observed, however, if the electron beams are properly 
polarized. 

Double chargino production e^e^ — > X X~ in. supersymmetry has been studied 
in"^^*. The diagram is basically the same as (a) and (b) in Fig. 37, with the W 
replaced by x ^^nd the neutrino by a sneutrino. Recall that lepton number violation 
in the sneutrino sector implies Majorana neutrinos^^'"' . It turns out that the limits 
from render double chargino production cross sections too small. Other works 
on lepton number violating e~e~ collisions within supersymmetry can be found 
in3i5. 

As mentioned before, lepton number violating sneutrino mass terms, the am- 
plitude of QvPP and Majorana neutrino mass terms imply each other: if one of the 
three is present, the other two are there as well^^^'^^^. This leads to a splitting in 
the v-v system and therefore to lepton number violating sneutrino-anti-sneutrino 
mixing'^^^, whose parameters depend heavily on SUSY parameters, and whose ob- 
servation is usually a very challenging task'^^''. 



8. Summary 

Neutrino-less double beta decay experiments are much more than neutrino mass 
experiments, their importance is much broader and deeper. The violation of lepton 
and baryon number is a rather generic feature of theories beyond the Standard 
Model, and searches for Ov/S/S or proton decay are probes of fundamental physics 
related to high energies, with a variety of important consequences in particle physics 
and cosmology. The significance of the decay is underlined by the excessive list of 
references provided in this review. 

In the next 20-30 years Oi//3/3 will be the only realistic probe to test the conser- 
vation of lepton number. The existing and upcoming results, when interpreted in 
terms of a specific particle physics scenario, allow to constrain a variety of impor- 
tant parameters, some of which can only be probed by Ovfi/S, others can also be 
tested in different and complementary experiments. Best motivated is presumably 
the standard interpretation of light neutrino exchange, where the inverted mass 
ordering will begin to be tested within this decade. Quasi-degenerate neutrinos will 
generate a signal, and should in this case be detectable also in direct searches and 
cosmological observations. This would be the ideal case to identify the mechanism. 
There are however many non-standard interpretations of 01^/3/3, the most frequently 
discussed mechanisms for the decay arc summarized in Table 14. The unambigu- 
ous determination of the underlying mechanism is in general less straightforward 
than for quasi-degenerate neutrinos. However, naive estimates show that an effec- 
tive mass of order 0.1 eV is associated with an amplitude that corresponds to the 
amplitude for exchange of TeV scale heavy particles. This energy scale has a vari- 
ety of potential effects in currently running particle physics experiments, such as 
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Table 14. Most important mechanisms for neutrino-less double beta decay. Given are the absolute value of the 
amplitude {q ~ 0.1 GeV is the momentum exchange for long-range processes) with the particle physics parameter 
written in bold face. The current limits on these quantities are provided, and tests to identify the mechanism by 
other means are indicated. RHC denotes right-handed currents, fi, stands for i?-parity violation and there are several 
Majoron variants. 



mechanism 



amplitude and 
particle physics parameter 



current limit 



test 



light neutrino exchange 



0.5 eV 



oscillations, 
cosmology, 
neutrino mass 



heavy neutrino exchange 




si. 

Mi 




2 X 10"* GeV-i 


LFV, 
collider 


heavy neutrino and RHC 


F W 


V 


2 


4 X 10~i« GeV^s 


flavor, 




collider 



Higgs triplet and RHC 



G% 



10-15 GeV-1 



flavor, 
collider 
e~ distribution 



A-mechanism with RHC 



F q 



flavor, 

1.4 X 10-10 GeV-2 coUider, 

distribution 



Tj-mechanism with RHC 



6 X 10- 



flavor, 
collider, 
e~ distribution 



\^ — 



short-range fi 



7 X 10-1* GeV-5 



collider, 
flavor 



long-range fi. 



sin 20" Ai3,A;,3 



2 X 10-13 GeV-2 
1 X lO-" GeV-3 



flavor, 
collider 



Majorons 



<x|<3x>l or Kgx>| = 



10-*... 1 



spectrum, 
cosmology 



LHC, lepton flavor violation, FCNC, etc. It should be noted that, though some 
progress was made in recent years, high precision physics with Ovpp will presum- 
ably not be possible: nuclear matrix elements are unlikely to be known with more 
than 20% precision. Currently, one has to accept the (shrinking) 0(1) ranges of 
NME calculations and perform analyses of Oi^/3/3-rcsults keeping them in mind. 

An impressive number of upcoming experiments promises an exciting future 
for the field. Multi-isotope determination of 01^(3/3 with different experimental ap- 
proaches will be possible and is crucial in order to make an unambiguous claim 
of observation, help clarifying the nuclear matrix element calculations, and dis- 
tinguish the different mechanisms. In order to identify the origin of ncutrino-lcss 
double beta decay (the "inverse problem" of 0iy(i/3) three different possibilities exist: 
via effects in other observables, via exploring the decay products, and via nuclear 
physics effects. After the violation of lepton number is established, a highly inter- 
esting physics program of identifying the underlying mechanism and its origin will 
be possible. 
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